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Abstract 

We present an overview of recent efforts to calculate the interactions among hadrons using lattice 
QCD. After outlining the techniques that are used to extract scattering parameters, we detail the 
latest calculations of meson-meson scattering, baryon-baryon scattering and multi-meson systems 
obtained with domain-wall valence quarks on the staggered MILC lattices by the NPLQCD col- 
laboration. Estimates of the computational resources required to achieve precision results in the 
baryon sector are presented. 



* Review prepared for the International Journal of Modern Physics C. 
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I. INTRODUCTION 



Lattice Quantum Chromodynamics (QCD) is emerging from — what should be considered 
to be — a thirty year research and development (R and D) phase, into the production phase, 
where it will produce precise and accurate calculations of non-perturbative strong interaction 
quantities. The last five years has seen the field evolve from an era of quenched calculations 
in small spatial volumes at large lattice spacings, to the present day, where fully-dynamical 
calculations with chiral symmetry at finite lattice-spacing are standard, the lattice volumes 
have a spatial extent of L ~ 2.5 — 3.5 fm and lattice spacings are 0.13 fm. Unfortunately, 
it is still the case that the light-quark masses, ruq, are larger than those of nature, with 
typical pion masses ~ 300 MeV. However, the next five years will see calculations at 
the physical light-quark masses, ~ 140 MeV, in large volumes, 6 fm, and at small 
lattice spacings, 0.06 fm become standard. This impressive evolution will continue with 
the increasing computer power dedicated to Lattice QCD calculations, and development of 
new and faster algorithms with which to generate lattice configurations, and to compute 
quark propagators and hadronic observables. 

The impact on nuclear physics of being able to reliably compute strong interaction quan- 
tities directly from QCD, the underlying theory of the strong interactions, cannot be over- 
stated. While it is important to recover what is known experimentally to high precision, 
that is obviously not the main objective, and certainly not a good reason to pursue such 
an effort. The reason for investing (both money and careers) in this area is to be able to 
calculate quantities of importance that cannot be accessed experimentally, or which can 
be measured with only limited precision. Nuclear physics is an incredibly rich field, whose 
phenomenology has been explored for decades. However, there is still little understanding 
of the connection to QCD and the basic building blocks of nature, quarks and gluons. Such 
a connection will be firmly established with lattice QCD, and as a consequence, the calcu- 
lation of quantities that are not experimentally accessible will finally become possible. On 
a more academic level. Lattice QCD will allow for an exploration of how nuclei and nuclear 
interactions depend upon the fundamental parameters of nature, and an understanding of 
the fine-tunings that permeate nuclear physics will finally be translated into (the possible) 
fine-tunings of the light-quark masses. 

Two important examples of how Lattice QCD calculations can impact nuclear physics are 
in the evolution of a supernova and in the structure of nuclei. In the evolution of supernova, 
one of the key ingredients that determines whether a supernova evolves into a black-hole 
or a neutron star is the nuclear equation of state (NEOS). The NEOS is determined by 
the dynamical degrees of freedom and their interactions. The degrees of freedom at a given 
density are determined by the mass of the hadrons in the medium, which depends upon their 
interactions. For instance, the mass of the S~ is expected to be significantly less than in 
free-space due to its interactions with neutrons, however, the precise mass-shift is uncertain 
due to the model dependence of existing calculations [1]. As the S~ carries the same charge 
as the electron, it can become energetically favorable to have the nuclear material composed 
of neutrons, protons and S~'s, as opposed to just neutrons and protons, due to the location 
of the fermi-levels. Experimentally, little is known about the S~-neutron interaction, and 
a precise Lattice QCD calculation of this interaction, and others, will greatly reduce the 
theoretical uncertainty in this particular potential contribution to the NEOS in hadronic 
systems with densities of a few times that of nuclear matter. As regards the structure of 
nuclei, we now have refined many-body techniques, such as Greens function Monte-Carlo 
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(GFMC) [2] with which to calculate the ground states and excited states of light nuclei, 
with atomic number 14. Using only the modern nucleon-nucleon (NN) potentials that 
reproduce all scattering data below inelastic thresholds with x^/dof ~ 1, such as AVig [3], 
ones fails quite dramatically to recover the structure of light nuclei. The inclusion of a three- 
nucleon interaction greatly improves the predicted structure of nuclei [2] . Lattice QCD will 
be able to calculate the interactions of multiple nucleons, bound or unbound in the same way 
it can be used to determine the two-body scattering parameters. For instance, a calculation 
of the three-neutron interaction will be possible. 

A. Aspects of QCD and Nuclear Physics 

The structure and interactions of all nuclei are determined completely by QCD and by 
the electroweak interactions, primarily electromagnetism. The Lagrange density of QCD 
is written in terms of the quark and gluon fields, and is manifestly invariant under local 
SU{3)c transformations. At short-distances the coupling between the quarks and gluons, and 
between the gluons and themselves, a<j(Q), is small and allows for processes to be computed 
as an asymptotic series in a s{Q)- asymptotic freedom. However, as the typical length scale of 
the process grows, as{Q) becomes large, and perturbative calculations fail for 1 GeV^. 
At long-distances, and hence low energies, the appropriate degrees of freedom are not the 
quarks and gluons that the QCD Lagrange density is written in terms of, but the hadrons, 
such as pions and nucleons. It is the properties of collections of nucleons, hyperons and 
mesons that defines the field of nuclear physics, a field that, simply put, is the exploration 
of the non-perturbative regime of QCD. In this low-energy regime, the approximate chiral 
symmetry of the QCD Lagrange density is spontaneously broken by the vacuum, giving 
rise to pseudo-Goldstone bosons and the approximate isospin symmetry. With such rich 
dynamics present in what naively looks like a simple theory it is not surprising that relatively 
little progress has been made in determining the properties and interactions of nuclei, even 
the simplest nuclei, directly from QCD. Understanding and calculating the properties and 
interactions of nuclei is further complicated by the fact that nuclear physics, and hence QCD, 
exhibits one or more fine-tunings. The values of the light quark masses, the scale of chiral 
symmetry breaking, A^, and the electromagnetic coupling constant take values such that 
the scattering lengths in the two-nucleon systems are both unnaturally large compared with 
the range of the nuclear interaction. Further, the location of energy-levels in carbon and 
oxygen are such that the triple-a process can proceed to produce enough carbon for us to 
be writing this review. Of course, one looks toward an anthropic explanation (for a lengthy 
discussion of the anthropic principle see Ref. jSj) of these fine-tunings, but at present we 
have no understanding of how the fine-tunings exhibited in nuclear systems translate into 
fine-tunings in the fundamental parameters of nature (if in fact they are fine-tuned). 

B. Aspects of Lattice QCD 

The only known way to solve QCD in the low-energy regime is numerically using Lattice 
QCD. Lattice QCD is a technique in which Euclidean space correlation functions are calcu- 
lated by a Monte-Carlo evaluation of the Euclidean space path integral. The calculations 
are performed in Euclidean space so that contributions to any given correlation function 
that have a large action are exponentially suppressed. This is in contrast with Minkowski 
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space in which large action contributions result in a complex phase which will average to 
an exponentially small contribution with nearby paths. Space-time is discretized with the 
quarks residing on the lattice sites, and the gluon fields residing on the links between lattice 
sites. The lattice spacing, b, (the distance between adjacent lattice sites) is required to be 
much smaller than the scale of chiral symmetry breaking (or any physical length scale in the 
system) for its effects to be small. The effects of a finite lattice spacing can be systematically 
removed by combining lattice calculations of the same correlation functions at several lattice 
spacing with the low-energy effective field theory (EFT) constructed to explicitly include 
the lattice spacing. The EFTs describing such calculations are somewhat more complicated 
than their continuum partners as they must reproduce that matrix elements of the Symanzik 
action constructed with higher dimension operators induced by the lattice spacing[5]. While 
the action lacks Lorentz invariance and rotational symmetry, it does possess hypercubic sym- 
metry. As computers have finite memory and performance, the lattice volumes are finite in 
all four space-time directions. Generally, periodic boundary conditions (BC's) are imposed 
on the fields in the space-directions (a three torus), while anti-periodic BC's are imposed 
on the fields in the time-direction, which in many cases is much larger than the space- 
directions ^. For the calculations we will be discussing in this review, the lattice volumes are 
large compared with the Compton wavelength of the pion, and deviations in single particle 
properties from their infinite volume values are exponentially small, generically ~ g-"iiri 
Finally, the cost of performing lattice calculations increases dramatically with decreasing 
quark mass, and presently, calculations at the physical quark masses, ~ 140 MeV are 
not possible. Currently, lattice calculations are performed at unphysical values of the quark 
masses, and the light quark mass dependence of the observable of interest, which can be 
determined perturbatively in the low-energy EFTs, is used to extrapolate to the physical 
light quark masses. Therefore, the practical situation with current Lattice QCD calculations 
is that they are performed at finite lattice spacing, within finite volumes and at unphysical 
quark masses. The appropriate EFT (e.g. X"PT, heavy-baryon-x-PT) is then used to ex- 
trapolate to the infinite volume, continuum limit of QCD where the parameters of nature 
reside. 

C. Diversion into Sociology 

Many of the members of the NPLQCD collaboration (of which the authors are a part of 
and with the logo displayed in Fig. [T]) spent years developing EFTs for nuclear physics, 
and more generally performing calculations of strong interaction processes in the low-energy 
regime. Around 2005, it started to become apparent to us that in order to make progress in 
low-energy QCD, the approximate chiral symmetries of QCD were, in many cases, no longer 
sufficient to make predictions for quantities at the precision required to impact experiments, 
or to shed light on fundamental questions in strong interaction physics. Basically, the number 
of counterterms that were commonly appearing in calculations beyond leading order (LO) 
or next-to-leading order (NLO) were large, and in many instances exceeded the number of 
experimental measurements that were available or would become available in the near future. 
We decided that the only path forward in this area was to use Lattice QCD to calculate these 



A linear combination of propagators generated with periodic and anti-periodic boundary conditions can 
be used to effectively double the length of the time-direction 
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counterterms. The focus of the NPLQCD collaboration is to calculate nuclear reactions from 
Lattice QCD, and our initial efforts have focused on calculating the elastic scattering of the 
lightest hadrons. It is these calculations that are the focus of this review. 



We wish to stress that we would not have been able to perform the calculations in the 
time frame that we have if it were not for the Lattice QCD infrastructure that existed in 
2005 and continues to exist and grow in the United States, namely USQCD and the efforts of 
the lattice group at Jefferson Laboratory. The USQCD collaborative effort ^ has allowed for 
the large scale production, and more importantly for us, sharing of lattice resources. The 
MILC collaboration has produced, and continues to produce, multiple staggered (Kogut- 
Susskind) lattice ensembles at different quark masses, lattice spacing and lattice volumes. 
Further, NPLQCD has shared resources with the LHP collaboration ^. Another significant 
contribution to our program is the software package Chroma [6l [7], developed by Robert 
Edwards and his team at Jefferson Laboratory, which we have used extensively to perform 
our calculations. 

II. FORMAL ASPECTS OF LATTICE CALCULATIONS FOR NUCLEAR 
PHYSICS 

In this section we highlight the formal aspects of Lattice QCD that have direct bearing on 
the extraction of scattering and reaction processes in multi-hadron systems. Intricate details 
of Lattice QCD calculations are beyond the scope of this review. In later sections, however, 
we discuss specific aspects of Lattice QCD calculations that are relevant to the discussions 
at hand. For even finer details, we refer the reader to a number of excellent texts on the 
subject [U [9l [To]. Throughout, we will be discussing lattices with cubic symmetry in the 
spatial directions, each of length L in lattice units. The time direction is assumed to be 
much larger than the spatial directions so that all calculations are assumed to be at zero 
temperature, T = 0. Further, while the lattice spacing is conventionally denoted by "a", 
we will use "b" in order to avoid conflict with the notation for scattering lengths which are 
universally denoted by "a" . 



^ http://www.usqcd.org/" 

^ http://www.jiab. org~dgr/lhpc/ 




FIG. 1: The NPLQCD logo. 
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A. Physics From Euclidean Space Correlation Functions 

In the continuum limit, the Euchdean space correlations functions (suitably Fourier trans- 
formed) are the sums of exponential functions, whose arguments are the product of Euclidean 
time with the eigenvalues of the Hamiltonian in the finite- volume associated with eigenstates 
that couple to the hadronic sources and sinks. At large times, the correlation function be- 
comes a single exponential dictated by the ground state energy and the overlap of the source 
and sink with the ground state. As an example, consider the pion two-point function, (t) 
generated by a source (and sink) of the form 7r"'"(x, t) = m(x, )f:)75(i(x, t), 

a+(t) = <o| ^"(^'^) ^^(O'O) |o) = E ^^(O'O)) • (1) 

X X 

The sum over all lattice sites at each time-slice, t, projects onto the spatial momentum 
p = states. The source 7r"'"(x, t) not only produces single pion states, but also all possible 
states with the same quantum numbers as a single pion. More generally, the source and 
sink are smeared over lattice sites in the vicinity of (x, t) to increase the overlap onto the 
ground state and lowest-lying excited states. Translating the sink operator in time via 
7r+(x,t) = e^*7r+(x,0)e"^* gives ^ 

CMt) = E ^ E vr-(x,0)|n)(n|vr+(O,0)|0) ^ — . (2) 

At finite lattice spacing, the correlation functions for Wilson fermions remain the sum over 
exponential functions, but for particular choices of parameters used in domain-wall fermions, 
the correlation functions exhibit additional sinusoidally modulated exponential behavior at 
short-times (with a period set by the lattice spacing). 

It is straightforward to show that the lowest energy extracted from the correlation function 
in Eqs. ([T]) and ^ correspond to the mass of the Tr+ (and more generally the mass of the 
lightest hadronic state that couples to the source and sink) in the finite volume. Finite 
volume effects are exponentially suppressed [11], therefore, while every lattice calculation is 
performed in a finite volume, as long as the volume is large compared to the pion mass, the 
mass of stable single particle states can be extracted with high accuracy. 

Once any correlation function, e.g. C^+(t), is calculated, the most common objective is 
to extract the argument of the exponential function that persists at large times. One way to 
do this is to simply fit the function over a finite number of time-slices to a single exponential 
function. A second method that is perhaps a bit more "pleasing to the eye" , is to form the 
effective mass function, e.g. 

where both t and Meff.(t) are in lattice units. At large times, Meff.(t) becomes a constant 
equal to the mass of the lightest state contributing to the correlation function ^. 



^ We assume the absence of external electroweak fields that exert forces on hadrons in the lattice volume. 
^ This is obviously the most simplistic approach one can take to this problem. To extract the ground state 
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B. Hadronic Interactions, the Maiani-Testa Theorem and Liischer's Method 

Extracting hadronic interactions from Lattice QCD calculations is far more complicated than 
the determination of the spectrum of stable particles. This is encapsulated in the Maiani- 
Testa theorem [13], which states that S-matrix elements cannot be extracted from infinite- 
volume Euclidean-space Green functions except at kinematic thresholds ^. This is clearly 
problematic from the nuclear physics perspective, as a main motivation for pursuing Lattice 
QCD is to be able to compute nuclear reactions involving multiple nucleons. Of course, it 
is clear from the statement of this theorem how it can be evaded, one computes Euclidean- 
space correlation functions at finite volume to extract S-matrix elements, the formulation 
of which was known for decades in the context of non-relativistic quantum mechanics [H] 
and extended to quantum field theory by Liischer [TSl IIB] • Liischer showed that the energy 
of two particles in a finite volume depends in a calculable way upon their elastic scattering 
amplitude and their masses for energies below the inelastic threshold. As a concrete example 
consider tt^tt^ scattering. A tt^^tt"'" correlation function in the Ai representation of the cubic 
group [IT] (that projects onto the s-wave state in the continuum limit) is 

|p|=P x,y 

In relatively large lattice volumes the energy difference between the interacting and non- 
interacting two-meson states is a small fraction of the total energy, which is dominated by 
the masses of the mesons. In order to extract this energy difference the ratio of correlation 
functions, GT,+^+{p,t), can be formed, where 

G..MP,t) ^ f"-;;;^^''!, - E A e-^^" * , (5) 

and the arrow denotes the large-time behavior of G^+7r+. The energy eigenvalue. En, and 
its deviation from the sum of the rest masses of the particle, AEn, are related to the center- 
of-mass momentum p„ by 

AEn = En - 2m^ = 2^/ pI + ml - 2m^ . (6) 

To obtain pcot6{p), where 6{p) is the phase shift, the square of the center-of-mass momen- 
tum, p, is extracted from this energy shift and inserted into [HI [151 [13 [H] 

pcot*W^^s((g)') . (7) 

which is valid below the inelastic threshold. The regulated three-dimensional sum is [19] 

|j|<A 

S(x) ^ Y^J^ - 47rA , (8) 



and excited states one can implement the method of Liischer and Wolff [12] in which the correlation 
functions resulting from different sources and sinks are calculated. The resulting matrix of correlation 
functions is diagonalized, and the effective mass function for each resulting eigenvalue can be used to 
extract the spectrum. 

^ An infinite number of infinitely precise calculations would allow one to circumvent this theorem. 



9 



where the summation is over all triplets of integers j such that | j | < A and the limit A ^ oo 
is implicit. Therefore, by measuring the energy-shift, AEn, of the two particles in the 




FIG. 2: The function 8(77) vs. rj, defined in Eq. ([s]), has poles only for 77 > 0. 



finite lattice volume, the scattering phase-shift is determined at AEn- In the absence of 
interactions between the particles, |pcot5| = 00, and the energy eigenstates in the finite 
volume occur at momenta p = 2nj/L. It is important to re-emphasize that this relation 
in Eq. ([T]) is valid relativistically [151 US]- Perhaps most important for nuclear physics 
is that this expression is valid for large and even infinite scattering lengths [12]. The only 
restriction is that the lattice volume be much larger than the range of the interaction between 
the hadrons, which for two nucleons, is set by the mass of the pion. 

For the scattering of two nucleons, the scattering length is known to be unnaturally large 
at the physical pion mass, and therefore, the relation in Eq. ([t]) will need to be used to 
extract the scattering parameters. For systems that are not finely-tuned, such as the n~^7r~^ 
system, an expansion in the volume can be used. In the large volume limit {L ^ |a|) 
the energy of the two lowest-lying continuum states in the Ai representation of the cubic 
group [I?] are [T51 1 



AEq 
AEi 



Aira 
ML2 ~ 



1 - Ci— + C2 



+ 



+ 0{L- 



12 tan (5o 



[l + c[ tan 60 + c'2 tan^ 5o + ... ] + 0(1' 



(9) 



where 60 is the s-wave phase-shift evaluated at p = 27i/L. The coefficients, which result from 
sums over the allowed momenta |15l HE] in the finite cubic volume, are ci 
C2 = +6.375183, c[ -- 

physics sign convention for the scattering length.) In addition, for a > with an attractive 



-2.837297, 

-0.061367, and c'2 = —0.354156. (Note that we use the nuclear 
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interaction ^ a bound state exists with energy [19] (in the large volume limit) 

12 1 



A^_i = -4- 
M 



1 + ^ I e~^^ + 

-fL l-27(pcot5)' 



(11) 



where (pcot S)' = pcot 5 evaluated at = —7^. The quantity 7 is the solution of 



7 + pcot (5|p2=_^2 = , (12) 

which yields the bound-state binding energy in the infinite- volume limit. As expected, the 
finite volume corrections are exponentially suppressed by the binding momentum. This is 
consistent with the corrections to a single particle state where the lightest hadronic excitation 
is the zero-momentum two-particle continuum state, as opposed to a state containing an 
additional pion for, say, the finite volume corrections to the single nucleon mass ^. 

In the limit where L <C \{pcot5)~^\, which is a useful limit to consider when systems 
have unnaturally-large scattering lengths, the solution of Eq. ([T]) gives the energy of the 
lowest-lying state to be 

AE'o = J^J^ld-i + d2Lpcot5o + ... ] 

A^i = [ d[ + 4 ^pcot^o + - ] , (13) 

where the coefficients are di = -0.095901, ^2 = +0.0253716, d[ = +0.472895, 4 = 
+0.0790234 and where pcot^o in A_Eo is evaluated at an energy AE = -^j^ di, while 

pcot^o in AEi is evaluated at an energy AE = d[. The values of the df^ are deter- 
mined by zeros of the three-dimensional zeta-functions, and the expressions for AEi and 
AEi, excluding AE_i, are valid for both positive and negative scattering lengths. 

Recently, we have performed Lattice QCD calculations of the energy of multiple 7r"*"'s, 
from which both the two-body scattering parameters and three-body interaction were de- 
termined [22I |23]. In order to make use of the lattice calculations, the ground-state energy 
of n identical bosons in a finite volume was required. The calculation of the energy-shift of 
n identical bosons at 0{L~'^) is straightforward but tedious. The energy-shift of the ground 



^ The extension of the Chowla-Selberg formula to higher dimensions [20] gives 

S(-x2) ^ _27r2a; + 67re-2^=^ + ... , (10) 

for large x, where the ellipses denote terms exponentially suppressed by factors of e-2v^7ra:^ more. 
® The finite volume dependence of bound states has been explored numerically in Ref . [21J . 
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state is El 



AEoin,L) 




a x2 



TX LJ \Ti LJ 
+ (2n - 7)XJ + (5^2 - 41n + 63) K, 
X^ - eX^ J + (4 + n - ^2)^2 + 4(27 - 15n + n2)X /C 

+ (14n^ - 22772^ + 91972 - 1043)/: 



+ 



+ 



n 



n 



192 



ttL 



3(n 



(To + Ti + 



-3)X 



3/L6 
3) X 



+ 0(L-«) . 



(14) 



where the geometric constants that enter are 



X = -8.9136329 
C = 6.9458079 
Sms = -185.12506 



J 

To 



16.532316 
-4116.2338 



JC = 8.4019240 
Ti = 450.6392 



(15) 



and (fe) =n\ /{n — k)\/k\. The last term in the last bracket of Eq. (^14h is the leading relativistic 
contribution to the energy-shift. Deviations from the energy-shift of ra-bosons computed with 
non-relativistic quantum mechanics arise only for three or more particles as the two-particle 
energy-shift has the same form when computed in non-relativistic quantum mechanics and 
in quantum field theory [151 [16]. The renormalization-scale independent, but volume depen- 
dent, three-body quantity 



r/3(/x) + 3^3 -47r log (/xL) —S 

m \ / TT^m 



(16) 



was determined in recent Lattice QCD calculations [221 123], where rj3{n) is the coefficient 
of the three-body interaction in the Hamiltonian. In Eq. ( [l5| , iSms is the value of the 
scheme-dependent quantity S in the Minimal Subtraction (MS) scheme used to renormalize 
the theory. The naive dimensional analysis (NDA) estimate of rjg is rjg ~ l/{m.,^f^). The 
leading terms in Eq. ( 14 ) evaluated at n = 2 reproduce those terms shown in Eq. ^ . 



1. Is it a Bound-State or a Scattering-State ? 

Obviously it is important to understand what infinite-volume physics can be extracted from 
finite-volume calculations. One important question that arises is: if a negative energy shifted 
state is calculated on the lattice at finite-volume, does it correspond to a bound state or to a 
scattering state? Clearly, calculations of the same correlation function in multiple volumes 
will allow for the exponential volume dependence of a bound state ^ to be distinguished from 



^ Corrections to the bound-state pole-condition in terms of i cot 5 that depend exponentially upon the 
volume j2T) are equivalent to the corrections to the bound state mass [19] that have the same form. 
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the power-law volume dependence of a scattering state. However, one can make an educated 
guess about the nature of the state by the magnitude of the energy shift. Consider a simple 
system whose scattering amplitude is dominated at low-energies by the scattering length and 
effective range, as is the case for the scattering of two nucleons. The location of the states in 
the lattice volume are determined by the solution of Eq. ([T]). In Fig. [s] we show the graphical 




FIG. 3: The functions S{x) and pcot6 vs. x, where pcot6 is evaluated at p2 = ^x. The left 
panel is evaluated for scattering parameters a = +1.1 fm and r = 0.5 fm and all other scattering 
parameters vanish. The volume parameters are those of coarse MILC lattices with L = 20 (solid, 
blue), L = 40 (dashed, purple), L = 80 (dashed, tan), and L = 120 (dot-dashed, green) and 
b = 0.125 fm. The right panel is evaluated at the same parameters as the left except a = —1.1 fm. 
The intercepts of the two curves corresponds to the location of the energy-eigenstates in the finite 
volume. 

solution to Eq. ([T]) for two systems, one with a = +1.1 fm and r = 0.5 fm (left panel) and 
the other with a = —1.1 fm and r = 0.5 fm (right panel). One finds that for energies giving 

a value of x = (2^) less than xq = —0.0959006 {pcot6 < 0) the state wiU likely become a 
bound state in the infinite- volume limit. In contrast, states with energies giving a value of 
xq {pcot6 > 0) will likely become a continuum state in the infinite-volume limit. These 
statements are at best rules of thumb, and one can construct obvious exceptions. Further, 
one can imagine scattering parameters in p cot 6 that modify these rules. 

Numerical explorations of the spectrum of compact scalar QED containing a loose bound 
state have been performed in multiple volumes [21]. The results of that work confirm the 
large- volume behavior of the spectrum expected from Eq. ([7| (TH] and shown in Fig. [sj 

C. Finite-Lattice Spacing Artifacts and Mixed- Action Effective Field Theories 

As space-time has been discretized in the process of performing the Lattice QCD calculation 
of any given observable, there is a systematic deviation between the calculation and the 
actual physical value of the quantity. One anticipates that in the limit that the lattice spacing 
is much smaller than any physical length scale associated with the strong interactions, the 
finite-lattice spacing effects can be neglected for practical purposes. For a fixed lattice 
volume, this requires a lattice with a large number of lattice sites, and at present time 
typical lattice spacing are b ~ 0.125 fm on the coarse MILC lattices, b ~ 0.09 fm on the fine 
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MILC lattices, and b ~ 0.06 fm on the super-fine MILC lattices. It is clear that the lattice 
spacing of the coarse MILC lattices is not that much smaller than the scale of chiral symmetry 
breaking, ~ 1 GeV ~ 0.2 fm. Naively, such lattice spacing would clearly be problematic if 
the leading lattice-spacing artifacts entered at order 0{b) which is the case for naive Wilson 
fermions, but most calculations today are performed with improved-fermions, such as Clover, 
Domain- Wall, Staggered or Overlap fermions which all have lattice spacing artifacts entering 
at (9(6^) The work we will review is performed with a mixed-action. With the exception 
of the potentials between B-mesons, the staggered lattice discretization of the quark action 
is used to describe the sea-quarks, while the Domain- Wall lattice discretization is used for 
the valence quarks. The two discretizations have their advantages and disadvantages, and 
the mixed-action scheme has proved to be an effective way to perform the calculations that 
we are interested in. 

Once the lattice calculation has been performed and the quantities of interest, such as 
scattering lengths, have been extracted using the finite-volume techniques described above, 
the finite-lattice spacing effects must be removed. Ideally, one would calculate the same 
quantity at multiple lattice spacings, and then extrapolate to the continuum using the ap- 
propriate finite-lattice spacing expression [271 EH [29] computed in the low-energy EFT, such 
as mixed-action chiral perturbation theory (MAx-PT) At present, we do not have suffi- 
cient computer resources to perform calculations at multiple lattice spacings. We have just 
completed our work on the available coarse MILC lattices and we are currently calculat- 
ing with the fine MILC lattices. However, one of the very exciting results to emerge from 
MAx-PT is that when certain mesonic observables are written as an expansion in m^//^ 
(for two-flavors), and more generally Mm/ fu for three-flavors, as measured in the lattice 
calculation lattice spacing artifacts can be eliminated at low-orders in the expansion by 
a field redefinition, due to the good chiral symmetry properties of the domain-wall (more 
generally Ginsparg-Wilson) fermions used in the valence sector. The mixed-action correc- 
tions for the 7r+7r+ scattering length have been determined in Ref. [30]. It was demonstrated 
that when the extrapolation formulae for this system are expressed in terms of the lattice- 
physical parameters there are no lattice-spacing-dependent counterterms at 0{b'^), 0{b^) or 
0{m^b'^) ~ 0{b^). This was explained to be a general feature of the two-meson systems 
at this order, including the non-zero momentum states [HI]. There are additional lattice- 
spacing corrections due to the hairpin interactions present in mixed-action theories, but for 
domain-wall valence propagators calculated on the asqtad improved MILC gauge configu- 
rations, these contributions are completely calculable without additional counterterms at 
NLO. They depend only upon the valence meson and staggered taste-identity meson mass 
splitting [301 EI] which has been computed [32]. This allows for a precise determination 



There are chiral symmetry breaking corrections entering at 0{h), but they are exponentially suppressed, 
a measure of which is the value of the residual mass. 
^-'^ MAx-PT is the low energy effective field theory describing the dynamics of the pseudo-Goldstone bosons 
in the situation where the sea-quarks are governed by one lattice action, while the valence quarks are 
governed by a different one, in our case domain-wall valence quarks and staggered sea quarks. This is 
a partially-quenched theory that allows for the lattice spacing to be systematically included as a small 
expansion parameter, starting at 0{b^). 

We denote quantities that are computed directly from the correlation functions, such lattice- 
physical quantities. These are not extrapolated to the continuum, to infinite-volume or to the physical 
point . 



14 



of the mixed-action corrections to the scattering lengths at the various pion masses. In 
two-flavor MA^-PT (i.e. including finite lattice-spacing corrections) the chiral expansion of 
the scattering length at NLO is jSIj 



,/=2 



m: 



mz 



167r2/2 



3 log 



J" 
6ml 



where it is understood that and are the lattice-physical parameters [21] and 



A2 — ™2 



m: 



(17) 



(18) 



contains the leading lattice-spacing artifacts, u denotes a valence quark and j denotes a 



sea-quark, and isospin-symmetry is assumed in both the sea and valence sectors, 
is the mass of a meson composed of two sea (valence) quarks of mass mj 



m 



(muu) 



(m„) and the dots 

denote higher-order corrections to the meson masses. Clearly Eq. (|17|), which contains all 



0(m^6^) and C(6^) lattice artifacts, reduces to the continuum expression for the scattering 
length |33] in the QCD limit where A^^ ^ The three-flavor MA^-PT expression for 
K^K^ and K^t^^ scattering are somewhat more complicated, for obvious reasons and we 
refer the reader to Ref. [Mj . 

The fact that baryons do not have the same special status with respect to chiral symme- 
try as the pseudo-Goldstone bosons means that the field redefinitions that parametrically 
suppress finite-lattice spacing contributions to 7r+7r+ scattering, do not exist for baryons. 
Therefore one expects to see finite lattice spacing contributions to the scattering lengths 
for nucleon-nucleon, and hyperon-nucleon scattering at 0{lP') (and residual mass- type con- 
tributions). In order to achieve high precision calculations of baryon-baryon scattering, 
calculations on lattices with a small lattice spacing will be required. 



D. Potentials From Lattice Calculations 

One remarkable feature of nuclear physics is that one can understand and compute to reason- 
able accuracy the properties and interactions of nuclei working with an energy-independent 
two-nucleon potential alone. Phenomenologically, one finds that the three-nucleon interac- 
tion is required to improve agreement with experiment, as is the inclusion of meson exchange 
currents into electroweak matrix elements. But the fact remains that the two-nucleon po- 
tential is the dominant interaction in nuclei. 

There is a burning desire to construct a nucleon-nucleon potential directly from 
QCD, and hence from Lattice QCD. One can extract an energy-dependent and source/sink- 
dependent potential, defined at one energy (the energy of the two nucleons determined in 
the finite volume), however, this contains no more information than the phase-shift at the 
energy determined with Liischer's method. One such potential was calculated in quenched 
QCD in Ref [35] but, for the reasons mentioned above, it cannot be used as an input in 
nuclear calculations and it cannot be meaningfully compared to traditional nucleon-nucleon 
potentials. 



The counterterm is the same counterterm that appears in continuum x-PT. 

In this context, the word potential means an energy-independent potential. 



15 



A nucleon-nucleon potential may be defined from Lattice QCD calculations in the same 
way that potentials are determined from experimental measurements of the elastic scattering 
cross-section. A large number of calculations are performed, producing values for the phase- 
shift, along with an uncertainty, over a wide range of low-energies, and a potential is defined 
that minimizes the x^/dof in a global fit to the lattice calculations. At present, there is 
no practical program underway to perform such an analysis due to limited computational 
power. 



E. Statistical Errors Associated with Calculations of Observables in the Nucleon 
and Systems of Nucleons 

Lattice QCD calculations of quantities are performed by a Monte-Carlo evaluation of the 
path integral. As such, an important aspect of any such calculation is its statistical error 
and it is useful to consider a discussion presented by Lepage [36j. Consider an observable 
that is extracted from the correlation function {0{t)), such as a correlation function resulting 
from a pion source at time t = and a pion sink at time t. In any calculation there are a 
finite number of configurations on which to perform measurements and a finite number of 
lattice sites on each lattice. For statistically independent measurements of 0{t), the noise 
to signal ratio of this correlator behaves as I/a/ZV in the limit of large A^. However, the 
time-dependence of this ratio depends upon the specific correlation function that is being 
calculated. 

Consider the case where {0{t)) is the correlation function associated with n pions, each 
with a source of the form 7r"'"(x, t) = m(x, t)75(i(x, t), or one that is smeared over neighboring 
lattice sites, 

m)) = (^^^(O'O)) ) - Aoe—' , (19) 

where the interactions between the pions have been neglected in the large time behavior. 
The variance of this correlator is estimated to be 

X^vr-(x,t)j (^5^7r+(y,t)j (^7r+(0,0)j (^7r-(0,0)j ) - {9{t)f 

(^2-^2) , (20) 

where the large time behavior of the variance is dictated by the lightest intermediate state 
that can be formed from the propagators emerging from the sources associated with 6{ty6{t) 
and coupling to the sinks (annihilation diagrams are not included). It then follows that the 
ratio of noise to signal behaves as 



cr 



which is independent of time. (Here a/x denotes the ratio of standard deviation to mean.) 
That is to say that for the single pion correlator and correlators involving arbitrary numbers 
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of pions the errors are time- independent, and importantly do not exponentially grow with 
time. 

The situation is, unfortunately, not so pleasant for systems involving any number of 
baryons. For the case of a single proton, the correlation function has the form 

rW) = E {P\^,t)f (0,0)) A;,e-^^' , (22) 

X 

where an interpolating field that has non-vanishing overlap with the proton is ~ 
M°'^C75(i^M*''^eaf,c, where a, 6, c are color indices and z is a spin-index. The variance of this 
correlation function is 

X 

Ap2 e"^™-* - Ajo ^ Ap2 e"^*"-* , (23) 

and therefore the noise to signal ratio behaves as 



a (j{t) 1 



gl^mp-fm^jt _ ^24) 



X {Oit)) ^ 

More generally, for a system of A nucleons, the noise to signal ratio behaves as 

^ gA(mp-|m,)t ^25) 

X Vat 

Therefore, in addition to the signal itself falling as G ~ ^-Ampt^ ^-^^ noise associated with 



the correlator grows exponentially as in Eq. (25). 



III. OUR CURRENT TECHNIQUES AND RESOURCES 

The Lattice formulation of QCD is the perfect tool for evaluating the correlation functions 
required to extract physical observables, such as hadron masses and phase shifts. It both 
provides an ultraviolet regulator of the continuum field theory and converts functional inte- 
grals into regular integrals of very high dimension. In the continuum, the QCD path integral 
is 

Z = j VA^VijV^jj e/'^'K-i^M^^"'^''-'^[^M7M+-]^ + Ca.^.) (26) 

where A^ is the gauge field representing the gluons, F^^, is the gauge field strength and -0, ip 
are the fermion fields representing the quarks. is the covariant derivative which ensures 
gauge invariance and 7^ are matrices satisfying the Clifford algebra. The physical quantities 
in this theory can be calculated from correlation functions of operators O that are functions 
of the quantum fields (quarks and gluons). 

{0) = ^J VA^V^ViP O e/'^'-(-i^M^^'^"''-^[^M7M+™]^ + Ca.F.) (27) 

We can now discretize the continuum path integral introducing a discrete space time. In 
order to preserve gauge invariance the gauge fields are discretized as special unitary matrices. 
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FIG. 4: A two dimensional slice of the four dimensional lattice. 

SU (3), living on the links of the lattice (see Figure |4]). The discrete gauge action is given by 
the sum over all plaquettes Pfj,u{x) which are the product of the links U going around the 
elementary plaquettes of the lattice. 

S,iU) = /5 E - ^ReTrP^^(a;)') (28) 

with 

P^.u = U,{x)U,{x + fi)Ul{x + u)Ul{x) , (29) 

and P is the lattice gauge coupling. Taking the lattice spacing to zero, the above action 
becomes the continuum gauge action —J d^x\ (-^^j,(a;)) . This is the well known Wilson 
gauge |37] action. This discretization is not unique but it is the simplest. One can modify 
this discrete action by adding larger loops with coefficients appropriately chosen in order 
to achieve better convergence to the continuum limit, which is the ultimate goal of the 
calculation. 

The fermions, which live on the vertices of the lattice, present a more challenging problem. 
Naive discretization results in the so called fermion doubling problem, i.e. lattice fermions 
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come naturally in sixteen copies, too many for describing real QCD which has three light 
quarks (up, down and strange) and three heavy quarks (top, bottom and charm). The dou- 
blers can be avoided by several ingenious formulations of lattice fermions. Wilson fermions, 
which were introduced first eliminate the doublers by adding irrelevant dimension five 
operators in the action that lift the masses of the doublers, leaving only the light fermion 
in the spectrum. The price to pay is breaking of chiral symmetry and the introduction of 
lattice artifacts that scale as 0{b). Kogut-Susskind fermions [38] provide another way to 
remove some of the doublers and re-interpret the remaining four as four degenerate flavors. 
In this approach a U{1) chiral symmetry still remains unbroken and lattice artifacts scale as 
0{b^). Kogut-Susskind fermions become problematic when the required number of flavors 
is not a multiple of four (as is the case for QCD). In addition, the broken flavor and chiral 
symmetries introduce large lattice artifacts, although they scale as 0(6^). Finally, the so 
called domain wall fermions [23 SHI SI] and overlap fermions US] are fermionic actions 
that both preserve chiral symmetry at finite lattice spacing and are doubler free. Unfortu- 
nately, such formulations are computationally significantly more expensive. In all cases the 
lattice fermion action is of the form 



Sf = i)D{U)i) (30) 

where tp is the fermion "vector" and D{U) is a sparse matrix acting on this vector, that 
depends on the gauge field U. 

The partition function in the case of two quark flavors is 



= [HdU^ix) det{D{UyD{U)) e'^^^^K (31) 

^,x 

The integration over the quark fields, which are represented by Grassman numbers, can 
be done exactly. In addition, the quark matrix D{U) represents one flavor, however since 
detD{Uy = detD(U), the determinant det [D(Uy D{U)) represents two flavors. In the case 
of correlation functions, integrating out the quarks gives the following expression 

= ^/ n^f^M(^) (^(^yU) det {D{U)^D{U)) e-'^^^K (32) 

resulting in operators O that depend on the inverse of the quark matrix. The above manipu- 
lation is only valid in the case of two flavors of quarks (the up and the down) which both have 
the same mass, which is a good approximation to the low energy physics of QCD. A strange 
quark can be easily added by including det [D{UyD{U)Y^'^ in the partition function. 



The computation of Eq. (32 ) is the main numerical task faced in Lattice QCD calculations. 
The integral in Eq. ( 32 ) over the gauge fields is of extremely large dimensionality. Consider- 
ing that we have discretized QCD, which has a fundamental scale of ~ 1 fm (10~^^ cm), we 
need to work with a lattice that has a physical size much larger than 1 fm in order to control 



In certain cases such as the overlap fermions the matrix is not sparse but has sparse hke properties i.e. 
the matrix vector multiphcation is a "cheap" operation. 
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finite volume effects, and a lattice spacing much smaller than 1 fm in order to control the 
continuum limit. With moderate choices for the volume and the lattice spacing, a typical 
lattice size of 32^ is arrived at. Counting the color, flavor and spin degrees of freedom, the 
calculation involves ~ 10* degrees of freedom. The only way this computation can be done 
is by using Monte Carlo integration. Fortunately, the combination of the quark determinant 
and the gauge action, 

V{U) = ^det {D{U)^D{U)) e-^«(^) , (33) 

is a positive definite quantity which can be interpreted as a probability and hence importance 
sampling can be employed. The basic algorithm is to produce gauge field configurations 
{f/} with probability distribution P(f/) and then evaluate 

1 ^ 1 

i=l ^ ' 

At finite A^, the estimate of O is approximate, with an error that converges to zero as 



0(l/v A^). Both for the gauge field configuration generation and the evaluation of Eq. (34), 
the linear system of equations 

D\U)\m\D{U)\m\x = <\>, (35) 

needs to be solved where the dependence of the quark matrix on the quark mass m is made 
explicit. Since the quark matrix is sparse, iterative solvers such as conjugate gradient can 
be used. The condition number of the quark matrix is inversely proportional to the quark 
mass. Since the physical quark masses for the up and down quarks are quite small, the 
quark matrix has a large condition number. With current computer resources this linear 
system cannot be solved exactly at the physical quark mass. For that reason the calculation 
is performed at heavier quark masses and then extrapolated to the physical point. The vast 
majority of the computer time used in these calculations is devoted to the solution of this 
linear system both in the context of gauge field generation and in the later stage of the 



calculation of physical observables through Eq. (34). 

Realistic lattice calculations require quark masses that result in pion masses below 400 
MeV, allowing chiral effective field theories to be used with some reliability. In addition, a 
dynamical strange quark is required in order to guarantee that the low energy constants of 
the EFT match those of the physical theory. Although this task seems formidable, in the 
last several years there have been developments that make phenomenologically interesting 
calculations now possible. 



A. Domain Wall and Staggered Fermions 

The emergence of fermions that respect chiral symmetry |39l HQl HH HSl HS] on the lattice 
was one of the major recent developments in Lattice QCD. These formulations of lattice 
fermions allow us to reduce the lattice spacing errors and approach the continuum limit in 
a smoother manner. However, the cost of calculating with these fermions is an order of 
magnitude larger than any other variant of lattice fermions. In addition, the development of 
improved Kogut-Susskind fermion actions [lllll5] that significantly reduce the 0{b^) errors, 
allowed for cheap inclusion of quark loop effects in the QCD correlation functions computed 
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on the lattice. With this formulation, volumes with spatial extent as large as L ~ 3.5 fm 
are possible with light-quark masses as low as 1/lOth of the strange quark mass depending 
on available computing resources. However the fact that Kogut-Susskind fermions represent 
four flavors of quarks complicates calculations when two or one flavors are needed. From 
the operational point of view the problem is solved by introducing into the path integral 
the Kogut-Susskind determinant raised to the ^^//4 power (rooted), where rif is the desired 
number of flavors. The non-integer power of the quark determinant introduces non-localities 
in the lattice action. It has been argued, however, that the long distance physics that survives 
the continuum limit is not affected by such non-localities [IHl HTj HHl HQ], [501 |5T]. In addition, 
at finite lattice spacing, the pathologies arising in the Kogut-Susskind fermion formulation 
can be dealt with in staggered x-PT |371 HHl [521 [531 Ell [SH [56] . Although no rigorous proof 
exists, empirical evidence indicates that Kogut-Susskind fermions do describe the correct 
physics as long as the continuum limit is taken before the chiral limit [ST]. It should be 
noted that there are some members of the lattice community who believe that the rooted- 
staggered action is fundamentally flawed and its continuum limit does not correspond to 
QCD (for a summary of these arguments, see Ref. [57]). We disagree with these arguments, 
however we acknowledge that there is no proof that the continuum limit of the rooted-Kogut- 
Susskind action corresponds to QCD. All of the work that we present in this review based 
upon mixed-action calculations on the MILC lattice ensembles assumes that the continuum 
limit is, in fact, QCD. 

In our calculations we use Kogut-Susskind fermions to represent the QCD vacuum po- 
larization effects associated with the two light flavors (up/down quarks) and the somewhat 
heavier strange quark. This is done by using gauge configurations generated with the appro- 
priate Kogut-Susskind fermion determinants incorporated into the probability distribution 
that enters the path integral. Since this part of the computation is completely discon- 
nected from the calculation of correlation functions, we can use gauge fields generated by 
other collaborations. In our case we use the gauge configurations generated by the MILC 
collaboration |58j . 

For all external quarks we use domain wall fermions. Because of the chiral symmetry 
that domain wall fermions satisfy, all our correlation functions satisfy chiral Ward identi- 
ties, ensuring that the leading order chiral behavior is continuum-like. The small corrections 
appearing due to Kogut-Susskind fermions in the vacuum loops can be taken care of system- 
atically in X"PT [301 EH [3l] . Compared to calculations with Kogut-Susskind fermions in the 
valence sector, this formulation results in better control of the chiral behavior and possibly 
smaller discretization errors. This approach was first introduced by the LHP collaboration 
for the study of nucleon structure [591 [SOI lEIl [621 [63] . 

B. The Lattice Actions 

The gauge configurations used in our work were generated by the MILC collaboration using 
the one loop tadpole improved gauge action [61] where both (9(6^) and 0{g%^) errors are 
removed. This action is defined as 



where R^]x,iiy and C[U]x,ijluct denote the real part of the trace of the ordered product of 
SU(3) gauge links along 1x2 rectangles in the /x, v plane and the /x, a, — /x, —v, —a paths, 




(36) 
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respectively. The coefficients cq, ci, and C2 are determined in one loop tadpole improved 
perturbation theory [6l], and [3 = 10 /qq where go is the bare gauge coupling. For the 
fermions in the vacuum the Asqtad improved Kogut-Susskind action [HI H^l ESI ESI EZl EH] 
is used. This action is the Naik action jSH] (0(6^) improved Kogut-Susskind action), with 
smeared links for the one link terms so that couplings to gluons with any of their momentum 
components equal to n/b are set to zero. The form of the action is the following: 



S 



Asqtad 
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Yl {yfiix)^y,x+fL - v/i{x - fi)Sy^x+fi) x{y)- 

. x,y fi=l,A 



x,y 11=1,4 ^ 



+ 



(37) 



where the Vn and links are 
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+ Ul{x - 0)Ul{x - 2u)U^{x - 2P)U^{x - 2u + fi)U^{x - + fi) 

Spup{x) = Uu{x)Spp{x + u)Ul{x + fi) + 

+ Ul{x - i')Spp{x - u)Uy{x - z> + /t) 

Spupa{x) = Uy{x)S^,p„{x + V)UI{X + fl) + 
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(38) 

(39) 
(40) 



(41) 
(42) 
(43) 
(44) 



In all of the above mq is the tadpole coefficient which was determined self consistently by 
MILC, rfpi^x) = (— l)Ep<M^p and x the four integers describing the lattice coordinates. 
For the valence sector we use the five dimensional Shamir domain wall fermion action [lOl 
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with D^{x,x') the regular four dimensional Wilson fermion action, 



= (4 + M5)5,,x' - 5^ 
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(45) 



(46) 



and Lg is the extent of the 5th dimension. Hypercubic-smeared (HYP-smeared) [TOl [71], [721 
gauge links were used in Eq. ( 45 ) and ( 46 ) to improve chiral symmetry. The physical four 



dimensional quark fields appear as boundary modes at the surface of the five dimensional 
space when M5 lies in the interval (—2,0). The physical quark fields (g(x) and q{x) are 
related to the underlying 5D fermions by 



q{x) 
q{x) 
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(47) 



The parameter m in Eq. (45) is related to the physical quark mass as it introduces in the 
effective action a mqq term. Domain wall fermions in the infinite Lg limit poses an exact 
chiral symmetry when m vanishes. This symmetry transformation is 



^(x,s) ^ ^(x,s)e-^^^(")^(^) 



(48) 
(49) 



where T^(s) 
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' Ls-l 
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However, at finite Lg this chiral symmetry is explicitly broken by the coupling of left 
handed and right handed modes in the middle of the 5th dimension. As a result one can 
construct the following partially conserved axial vector current 



L,-l 



= - X] ^5{s)j^{x,s) 



(50) 



s=0 



where is the four dimensional conserved vector current that corresponds to the 4D Wilson 
fermion action. This current satisfies a Ward-Takahashi identity which in the flavor non- 
singlet case takes the form [H]: 



(51) 



A,{A''^{x)0{y)) = 2m{q{x)T^j,q{x)0{y)) + 

+ 2{q^p{x)T'^^^q^p{x)0{y)) + i{5-0{y)) 



23 



where 



= i{x. y - 1)^^ + mx, y)i±^ (52) 
are four dimensional fields constructed at the midpoint of the the 5th dimension. The 



Ward-Takahashi identity of Eq. (51) is the same as the continuum counterpart with just 
an additional term 2{qmp{x)T"''y5qmp{x)0{y)) . This term is there only at finite Lg and 
it is a measure of the explicit chiral symmetry breaking. At long distances this term is 
proportional to 2(g(x)r°75g(a;)0(?/)). Using the pseudo-scalar density as a probe operator 
0{y) the residual mass is defined as 

> (gmp(t)r"75gmp(t)g(0)r"75g(0)) 



^ ^ 1 y^^ (gmp(t)r°75g^p(t)g(0)r"75g(0)) 



to 



where to,tmax is the time interval where only the ground state pion contributes to the two 
correlators in the ratio. 



C. Mixed Action Tuning 

Because the valence and sea quark actions are different, the calculation is inherently partially 
quenched. In other words, the calculation violates unitarity. Unlike conventional partially 
quenched calculations, which become unitary when the valence quark mass is tuned to the 
sea quark mass, unitarity cannot be restored by tuning the valence quark mass. The next 
best option is to tune the valence quark mass in such a way that the resulting pions have the 
same mass as those made of the sea Kogut-Susskind fermions. In this case unitarity should 
be restored in the continuum limit, where the Uf = 2 staggered action has an SU{8)l <8) 
SU{8)r ® U{l)v chiral symmetry due to the four-fold taste degeneracy of each flavor, and 
each pion has 15 degenerate additional partners. At finite lattice spacing this symmetry 
is broken and the taste multiplets are no longer degenerate, but have splittings that are 
C(o;^6^) im HHl EHl ESI EH]- The domain wall fermion mass is tuned to give valence pions 
that match the Goldstone Kogut-Susskind pion This choice gives pions that are as 
light as possible, resulting in better convergence of the X"PT needed to extrapolate the 
lattice results to the physical quark mass point. This tuning was also done by LHPC 
collaboration [59l EQl EH ESI [H EH] . 



D. Method of Contractions 



Throughout our work, gauge invariant Gaussian smeared quark propagators centered around 
a single point were used. In order to facilitate the complicated Wick contractions of many 



For the flavor singlet current this term survives the infinite hmit and gives rise to the anomaly. 
This is the only Goldstone boson that becomes massless in the chiral limit at finite lattice spacing. 
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body interpolating fields, the following strategy was adopted. Noticing that the permu- 
tations needed for Wick contractions are in effect "scalar" all the contractions at the 
annihilation operator point (sink) were performed and all the color and spin indices were 
left open at the creation operator point (source). The resulting data are Fourier transformed 
(space indices at the annihilation operator point) and saved on disk. All the two-body (and 
N-body) correlation functions can then be constructed by appropriate contractions of the 
source spin and color indices on scalar machines (such as simple work stations). 

Code was constructed to automatically perform all permutations, keeping track of the 
signs associated with fermion exchanges, allowing the construction of complicated diagrams 
in a relatively simple and efficient manner. This approach works well when quark annihi- 
lation diagrams are absent. For this reason such processes (such as the 1=0 vrvr channel), 
have not been explored in our work to date. Another reason for avoiding such processes is 
that the mixed action approach has large artifacts in these channels. In the future we plan 
to address such processes using the so called all-to-all propagator approach [75] . 



E. Data Analysis 

Since Monte-Carlo integration is used to compute the relevant correlation functions, the 
statistical uncertainty must be carefully determined. The main observables extracted in all 
calculations presented in this review are energy levels and energy level differences. These 
results contain information about phase shifts, scattering lengths and the three body inter- 
action as discussed above. The extraction of energy levels is done by fitting the relevant 
correlation functions to a sum of exponentials (or appropriately hyperbolic cosine functions 
when anti-periodic boundary conditions in time are imposed). We performed correlated 
minimization fits that take into account the time correlations in the lattice data. In partic- 
ular, the relevant parameters, such as the masses and the amplitude each state contributes 
to the correlation function, are determined as those that minimize 

X\A) = [G{U) - F{U, A)] [Git,) - Fit,, A)] (54) 

where G{t) are the lattice two point correlation functions, F{t, A) are the fitting functions 
used, with A denoting the set of fitting parameters over which x^iA) is minimized, and 
Cij is the covariance matrix. The lattice two point correlation functions are determined as 
averages over Monte-Carlo samples Gfe(t): 

1 ^ 

Git) = j^Y.^kit) (55) 

k=l 

and 

1 ^ 

C^, = E [^'^(^^) - ^(^^)] [^'^(^^■) - ^(^^■)] ■ (56) 

The (standard) errors on the fitted parameters are determined by the boundaries of the error 
ellipsoid, which is defined by the locus of points where = Xmin + 1- (For a pedagogical 
presentation of fitting see the TASI lectures by D. Toussaint [77].) 



In parallel computing language "scalar" refers to operations that are not data parallel. 
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In computing scattering lengths, the procedure for determining the statistical errors is a 
little more involved due to the highly non-linear relation between the scattering length and 
the energy levels of the two particle system. First one is interested in the energy differences 
between the energy levels of the two particle system and the sum of the masses of the two 
free particles (similarly for the case of more than two particles). These energy differences 
can be determined in two ways. The simplest is the one described in the previous sections 
where ratios of correlators are constructed in such way so that these ratios are a sum over 
exponentials parametrized by the desired energy splittings (see for example Eq. ([6])). In this 
case Jackknife is used to determine the covariance matrix and then a correlated x^"fit is 
performed as described above. For a single elimination Jackknife, the covariance matrix of 
a ratio of correlators is 

N-1 ^ 

k=l 

where Rk is the desired ratio computed with the kth sample omitted from the full ensemble 
and R is the ratio computed on the full ensemble. 

Fitting correlation functions to extract the ground state energy requires fitting ranges 
that start at time separations from the source that are large enough so that excited states 
have negligible contributions. The determination of the minimum time separation that can 
be included in the fit is sometimes subjective. Hence a systematic error due to the choice of 
the minimum time separation in the fit is included. This error is determined by observing 
the variation of the extracted results as a function of the choice of fitting interval. The final 
errors include both systematic and statistical errors combined in quadrature. 

One way to reduce the systematic errors due to coupling to excited states is to use 
appropriate interpolating fields that reduce these couplings. This can be achieved by using 
smeared quark sources and sinks. We have used gauge invariant Gaussian smeared sources 
for the quark propagators which were implemented as 

X(x) = (l + ^V^) S{xo) (58) 

where 6{x) is a 3D delta function, is the covariant 3D Laplacian, and a; is a parameter 
controlling the width of the smearing. In the limit of N oo this becomes an exponential 
of the Laplace operator. While the value of uj was chosen to optimize the coupling to the 
nucleon, it also worked very well for the mesons. In addition, we used both point and 
smeared sinks, both of which produced comparable results. 

All of our lattice calculations are performed away from the physical pion mass. For that 
reason we have to rely on x-PT to extrapolate to the physical point. Since results from 
ensembles with different sea quark masses are un-correlated (i.e. the covariance matrix is 
diagonal), we performed un-correlated fits of our data to the x-PT formulas. In all cases 
we estimate a systematic error due to these fits and extrapolation. Finally, except for 7r"'"7r"'" 
and K^K^ scattering lengths where the leading order lattice spacing effects are calculable, 
lattice spacing errors are estimated by dimensional analysis. 

In order to avoid systematic errors due to scale setting, dimensionless quantities (appro- 
priate products or ratios of dimensional quantities) are used wherever possible. In the few 
cases we need to quote results in MeV we use the scale determined by the MILC collaboration 
and confirmed by our own determination using the pion decay constant . 



26 



An essential parameter for all of our calculations using domain wall fermions in the valence 
sector is the residual mass uires that measures the degree of explicit chiral breaking. The 
results of the residual mass calculation are presented in Table [1} 



TABLE I: The parameters of the MILC gauge configurations and domain-wah propagators used in 
this work. The subscript / denotes hght quark (up and down), and s denotes the strange quark. The 
superscript dwf denotes the bare-quark mass for the domain-waU fermion propagator calculation. 

Ensemble bmi buis brnf^-^ brris^^ 10^ x brrires 

2064f21b676m007m050 0.007 0.050 0.0081 0.081 1.604 
2064f21b676m010m050 0.010 0.050 0.0138 0.081 1.552 
206412 Ib679m020m050 0.020 0.050 0.0313 0.081 1.239 
2064f21b681m030m050 0.030 0.050 0.0478 0.081 0.982 
2896f2b709m0062m031 0.0062 0.031 0.0080 0.0423 0.380 



F. Resources 



The computational resources needed for the calculations presented here were obtained from 
several sources including the USQCD clusters at JLab and FNAL (pentium clusters with 
infiniband interconnect), Tungsten at NCSA (pentium infiniband cluster) and Mare-Nostrum 
in Barcelona, Spain (Power PC 970MP IBM with Myrinet interconnect) The total 
computer power that went into our results reviewed here is about 1.2Tfiop- years, an amount 
of time that is rather small by current standards. This amount of time was sufficient to 
achieve precision results in the meson sector. 

The results presented here are snapshots of an ongoing effort. Not all of them were 
obtained with the same statistics on a given ensemble. For details on exactly what went into 
these calculations the reader should refer to the original publications. We typically computed 
several propagators per configuration by shifting the location of the source, and have cases 
with as many as 24 propagators per configuration. We carefully monitored the variance 
reduction as more propagators were added, and concluded that correlation functions were 
almost statistically independent. The residual correlations were taken care of by averaging 
correlation functions on a given configuration (blocking) and proceeding with the statistical 
analysis of data discussed in Section HIE Further blockings over multiple configurations 
did not change the results. 



IV. TWO-BODY PHYSICS 

A. TTTT Scattering 
1. Introduction 

Pion-pion scattering at low energies is the simplest and best-understood hadron-hadron 
scattering process. Its simplicity and tractability follow from the fact that the pions are 



We thank Andrew Pochinsky and Balint Joo for essential help in optimizing our code for this machine. 
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identified as the pseudo-Goldstone bosons associated with the spontaneous breaking of the 
approximate chiral symmetry of QCD. For this reason, the low-momentum interactions of 
pions are strongly constrained by the approximate chiral symmetries, more so than other 
hadrons. The scattering lengths for nn scattering in the s-wave are uniquely predicted at 
LO in x-PT [79]: 

m^4=° = 0.1588 ; mW^f = -0.04537 , (59) 

at the charged pion mass. Subleading orders in the chiral expansion of the tttt amplitude 
give rise to perturbatively-small deviations from the tree level, and contain both calculable 
non-analytic contributions and analytic terms with new coefficients that are not determined 
by chiral symmetry alone |33l EQl H] . In order to have predictive power at subleading orders, 
these coefficients must be obtained from experiment or computed with Lattice QCD. 

Recent experimental efforts have been made to compute the s-wave vrvr scattering lengths, 
ai=° and ai=^: E865 [821 [83] {Ke^ decays), CERN DIRAC [M] (pionium lifetime) and CERN 
NA48/2 [85] {K^ vr^vr^vr^). Unfortunately, these experiments do not provide stringent 
constraints on a^^^. However, a theoretical determination of s-wave vrvr scattering lengths 
which makes use of experimental data has reached a remarkable level of precision [HSl IH7] : 

m^4=° = 0.220 ±0.005 ; mW^=^ = -0.0444 ± 0.0010 . (60) 

These values result from the Roy equations [881 Ell |90] , which use dispersion theory to relate 
scattering data at high energies to the scattering amplitude near threshold. In a striking 
recent result, this technology has allowed a model-independent determination of the mass 
and width of the resonance with vacuum quantum numbers (the a meson) that appears 
in the vrvr scattering amplitude [91]. Several low-energy constants of one- loop X"PT are 
critical inputs to the Roy equation analysis. One can take the values of these low-energy 
constants computed with Lattice QCD by the MILC collaboration J3T, "92] as inputs to the 
Roy equations, and obtain results for the scattering lengths consistent with the analysis of 
Ref. [81]. 

The first lattice calculations of tttt scattering were performed in quenched QCD [931 EH 

E^EsiEziEHiEsiiinnittniiiiniiMiinaMM ra, 

and the first partially-quenched QCD calculation of tttt scattering (the scattering length and 
phase-shift) was carried through by the CP-PACS collaboration, who exploited the finite- 
volume strategy to study 1 = 2, s-wave scattering with two flavors [uf = 2) of improved 
Wilson fermions |114] . with pion masses in the range m^r — 0.5 — 1.1 GeV. The first fully- 
dynamical calculation of the J = 2 tttt scattering length with three fiavors (n/ = 2 + 1) 
of light quarks was performed by NPLQCD using domain-wall valence quarks on asqtad- 
improved staggered sea quarks at four pion masses in the range ~ 0.3 — 0.5 GeV at a 
single lattice spacing, b ~ 0.125 fm [78] . 



2. Data Analysis and Chiral and Continuum. Extrapolation 

It is convenient to present the results of the lattice calculation in "effective scattering length" 
plots, simple variants of the effective-mass plots discussed above. The effective energy split- 
ting is formed from the ratio of correlation functions given in Eq. (|5]), 



28 



which in the hmit of an infinite number of gauge configurations would become a constant 
at large times that is equal to the lowest energy of the interacting 7r"'"7r"'" system in the 
volume relative to the non-interacting two-pion energy. At each time-slice, Ai?7r+7r+(^) is 
inserted into Eq. ([t]) (or Eq. (|9])), to give a scattering length at each time slice, a^+^+{t). 
It is customary to consider the dimensionless quantity given by the pion mass times the 
scattering length, 077+ 7r+, where m7r(t) is the pion effective mass, in order to remove 
scale-setting uncertainties. For each of the MILC ensembles that were analyzed, the effective 
scattering lengths are shown in Fig. |5} and the values of the pion masses, decay constants 
and 7T7T energy-shifts can be found in Refs. [751 1115] . The mixed-action corrections for the 
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FIG. 5: The effective vr+vr"'" scattering length times the effective vr mass as a function of time-slice 
arising from smeared sinks. The sohd black lines and shaded regions are fits with l-cr statistical 
uncertainties; the dashed lines are estimates of the systematic uncertainty due to fitting. 



2 TiTT scattering length have been determined in Ref. [30] and have been discussed above 
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With domain- wall fermion masses tuned to match the staggered Goldstone pion [TU [75] . 
one finds (in lattice units) A|„ = fe^A^ = 0.0769(22) [32] on the coarse MILC lattices with 
b ~ 0.125 fm. For each ensemble, m.,r a^^^ was determined and the chiral extrapolation 



formula in Eq. (17) was used to extract a value of the counterterm /^^^(/x = fn)- The results 
of the two-fiavor extrapolation to the continuum are given in Ref. |115] . 

Fitting to lattice data at the lightest accessible values of the quark masses will op- 
timize the convergence of the chiral expansion. While there are only four different 
quark masses in the data set, with pion and kaon masses of approximately (771,^,^^^) ~ 
(290, 580), (350, 595), (490, 640) and (590, 675) MeV, fitting all four data sets and then "prun- 
ing" the heaviest data set and refitting provides a useful measure of the convergence of the 
chiral expansion. Hence, in "fit A", we fit the /^^^(/7r)'s extracted from all four lattice 
ensembles (m^ ~ 290,350,490, and 590 MeV) to a constant, while in "fit B", we fit the 
Z^=^(/^)'s from the lightest three lattice ensembles (m^ ~ 290,350, and 490 MeV). In "fit 
C", we fit the /^^^(/7r)'s from the lightest two lattice ensembles (m^ ~ 290 and 350 MeV). 
Results are given in Table [Tij Taking the range of parameters spanned by fits A-C one finds: 

TABLE II: Results of the fits in two-flavor Mixed-Action x-PT. The values of m,r a^^^ correspond 
to the extrapolated values at the physical point. The first uncertainty is statistical and the second 
is a comprehensive systematic uncertainty. 

FIT /^=2(^ = m^r a^^^ (extrapolated) X^/dof 
A 6.43 ±0.23 ±0.26 -0.043068 ±0.000076 ±0.000085 1.17 
B 5.97 ±0.29 ±0.42 -0.043218 ± 0.00009 ± 0.00014 0.965 
C 4.89 ±0.64 ±0.68 -0.04357 ± 0.00021 ± 0.00022 0.054 



lififi = U) = 5.4 ±1.4 

al=^ = -0.04341 ± 0.00046 . (62) 

In Fig. [6] we show the results of our calculation, along with the lowest mass nj = 2 point 
from CP-PACS (not included in the fit). The tree-level prediction and the results of the two- 
fiavor fit described in this section are also shown. The experimental point shown in Fig. [6] 
is not included in the fit and extrapolation. It is interesting that the lattice calculations 
indicate little deviation from the tree level X"PT curve. The significant deviation of the 
extrapolated scattering length from the tree-level result is largely a consequence of fitting 
to MA^-PT at one-loop level. 

An important check of the systematic uncertainties involved in the chiral extrapolation is 
to perform the same analysis using three-fiavor MA^-PT [301 l3T] as both the real world and 
our lattice calculation have three active light fiavors. In addition to the computations needed 
for the two-fiavor analysis, it is necessary to determine masses and decay constants for the 
kaon and the rj. The Gell-Mann-Okubo mass-relation among the mesons is used to determine 
the 7] mass, which is not computed in this lattice calculation due to the enormous computer 
resources (beyond what is available) required to compute the disconnected contributions. 
This procedure is consistent to this order in the chiral expansion. 

The chiral expansion of the 7r+7r+ scattering length in three-fiavor mixed-action X"PT 
as well as the numerical values for the various ensembles are given in Ref. [115] . For the 



three-fiavor analysis, the pruning analysis gives the results shown in Table III Taking the 
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m I f 

FIG. 6: m,r a-l^"^ vs. 171^^/ f.,^ (ovals) with statistical (dark bars) and systematic (light bars) uncer- 
tainties. Also shown are the experimental value from Ref. [83] (diamond) and the lowest quark 
mass result of the nj = 2 dynamical calculation of CP-PACS |114j (square). The blue band cor- 
responds to a weighted fit to the lightest three data points (fit B) using the one-loop MAx-PT 
formula in Eq. ( |17| ) (the shaded region corresponds only to the statistical uncertainty). The red 
line is the tree-level result. The experimental data is not used in the chiral extrapolation 

fits. 

TABLE III: Results of the NLO fits in three-flavor Mixed-Action x-PT. The values of ai='^ 
correspond to the extrapolated values at the physical point. The first uncertainty is statistical and 
the second is a comprehensive systematic uncertainty. 

FIT 32(47r)L^=2(/i = /^) a^^f (extrapolated) xVdof 

D 7.09 ±0.23 ±0.23 -0.042992 ± 0.000076 ± 0.000077 0.969 
E 6.69 ±0.29 ±0.39 -0.04312 ± 0.00009 ± 0.00013 0.803 
F 5.75 ±0.63 ±0.64 -0.04343 ± 0.00021 ± 0.00021 0.073 

range of parameters spanned by fits D-F one finds: 

32(47r)L^f(/i = U) = 6.2 ± 1.2 

a{=^ = -0.04330 ± 0.00042 . (63) 

3. Systematic Uncertainties 

There are many sources of systematic uncertainty to be quantified; there are lattice-spacing 
artifacts that arise at 0(m^6^); there are exponentially-suppressed finite- volume effects; 
there are effects due to residual chiral symmetry breaking; there are generic two-loop effects 
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due to the truncation of the chiral expansion at one loop; and finally, there are range 
corrections that enter at O {L~^) in Eq. (|9|. All of these effects have been taken into 
account |115j . It is noteworthy that the residual mass turns out to be one of the leading 
systematic errors in our analysis of 7r+7r+ scattering. 

It is worth emphasizing that the calculations have exact isospin symmetry, as do the 
extrapolation formulas used to analyze the results. The conventional discussion of the 
scattering length is in the unphysical theory with e = and m„ = rud = m, with 
m„ = = 139.57018 ± 0.00035 MeV and U = fn+ = 130.7 ± 0.14 ± 0.37 MeV. Hence 
m7r+//7r+ = 1.0679 ± 0.0032, where the statistical and systematic uncertainties have been 
combined in quadrature. The results of the lattice calculations are extrapolated to this value. 
The leading contribution to isospin breaking in vrvr scattering is due to the electromagnetic 
interaction, and this has been studied extensively Such contributions must be removed 
from the experimentally-determined scattering amplitude in order to make a comparison 
with the strong-interaction calculations. Isospin breaking due to the difference in mass of 
the light quarks occurs at next-to-leading order in the chiral expansion, and is expected to 
be small, as is its contribution to m'^, — m?n- 



4- Discussion 



J=2 



The prediction for the physical value of the J = 2 tttt scattering length is mj^a^^^ 
—0.04330 ± 0.00042, which agrees within uncertainties with the (non-lattice) determination 
of CGL [86]. In Table IV and Fig. PHwe offer a comparison between various determinations 



TABLE IV: A compilation of the various calculations and predictions for the 1 = 2 tttt scattering 
length. 



X-PT (Tree Level) -0.04438 

NPLQCD (2007) -0.04330 ± 0.00042 

E 865 (2003) -0.0454 ± 0.0031 

NPLQCD (2005) -0.0426 ± 0.0018 

MILC (2006)* -0.0432 ± 0.0006 

MILC (2004)* -0.0433 ± 0.0009 

CGL (2001) -0.0444 ±0.0010 



For discussions of the contributions of virtual photons to tttt scattering see Ref. |116l 11171 11181 1119] , and 
for very recent work on the electromagnetic contributions to the extraction of tttt scattering lengths from 
kaon decays see Ref. [120j and Ref. [121j 

The stars on the MILC results indicate that these are not lattice calculations of the / = 2 tttt scattering 
length but rather a hybrid prediction which uses MILC's determination of various low-energy constants 
together with the Roy equations), and the Roy equation determination of Ref. [86 (CGL (2001). 
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FIG. 7: Bar chart of the various determinations of the 1 = 2 vrvr scattering length tabulated in 
Table HVl See footnote EH 



B. Kvr Scattering 

1. Introduction 

In hadronic atoms, nature has provided a relatively clean environment in which to explore 
the low-energy interactions of charged hadrons. The electromagnetic interaction allows for 
oppositely-charged, long-lived hadrons to form Coulomb bound states. The locations of the 
energy-levels of these systems are perturbed by the strong interactions, while the lifetimes 
of the ground states are dictated by the strong interactions that couple the charged hadrons 
to lighter neutral ones. 

Study of the low-energy interactions between kaons and pions with -k^ bound-states 
allows for an explicit exploration of the three-flavor structure of low-energy hadronic in- 
teractions, an aspect that is not directly probed in vrvr scattering. Experiments have been 
proposed by the DIRAC collaboration [1225 study ttK atoms at CERN, J-PARC and GSI, 
the results of which would provide direct measurements or constraints on combinations of 
the scattering lengths. In the isospin limit, there are two isospin channels available to the 
TiK system, I = \ and I = \- The width of a t^'K^ atom depends upon the difference be- 
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tween scattering lengths in the two channels, F ~ (ai/2 — 03/2)^, (where ai/2 and 03/2 are the 
/ = ^ and / = I scattering lengths, respectively) while the shift of the ground-state depends 
upon a different combination, /\Eq ~ 2ai/2 + 03/2- Recently, the Roy-Steiner equations 
(analyticity, unitarity and crossing-symmetry) have been used to extrapolate high-energy 
T^K data down to threshold |123j . where it is found that 

(ai/2 - 03/2) = 0.269 ±0.015 , (ai/2 + 203/2) = 0.134 ±0.037 , (64) 

which can be decomposed to 7717^01/2 = 0.224 ± 0.022 and 771,^03/2 = —0.0448 ± 0.0077. (See 
also Refs. |124tll25] and Ref. |126j for a similar approach.) In addition, three-flavor X"PT has 
been used to predict these scattering lengths out to next-to-next-to-leading order (NNLO) 
in the chiral expansion. At NLO [ETl IT281IT29] . 

(ai/2 - 03/2) = 0.238 ± 0.002 , (ai/2 ± 203/2) = 0.097 ± 0.047 , (65) 

while at NNLO |130] mTrai/2 = 0.220 and m-,^a^/2 = —0.047 One must be cautious in 
assessing the uncertainties in these theoretical calculations, as one can only make estimates 
based on power-counting for the contribution of higher-order terms in the chiral expansion. 
There has been one determination of the t^^K^ scattering length in quenched QCD [131] , 
however, the chiral extrapolation of the scattering length did not include the non-analytic 
dependences on the light quark masses that are predicted by x-PT. 

As discussed above, recent work has identified in a model-independent way the lowest- 
lying resonance in QCD which appears in vrvr scattering Crucial to this development has 
been the accurate determination of the low-energy tttt scattering amplitude, including the 
recent Lattice QCD determination of the 1 = 2 scattering length [78] . A similar analysis has 
very recently been carried out for ttK scattering in the I = ^ s-wave in order to determine 
the lowest-lying strange resonance |132] . Improved accuracy in the low-energy nK scattering 
amplitude should be welcome to this endeavor. 

2. Analysis and Chiral Extrapolation 

It is useful to consider the dimensionless quantity of the reduced mass times the scattering 
length, /Ittk a^+i^+, where f^nKit), the "effective reduced mass" is constructed from the 
effective mass of the single particle correlators. For each of the MILC ensembles, the effective 
scattering lengths are shown in Fig. |8j while the results of the lattice calculation of the decay 
constants, meson masses, 7t~^K^ energy shifts and scattering lengths are given in Ref. |133] . 
The scattering lengths as a function of reduced mass are shown in Fig. |9} 

In SU(3) x-PT [531 11341 1135] at NLO, the expansion of the crossing even (a"^) and 
crossing odd (a~) scattering length times the reduced mass is known |127l 11281 1129] and 
reproduced in Ref. [133] . The counterterm L^;^(A) is a renormalization scale, A, dependent 
linear combination of the Gasser-Leutwyler counterterms 

L^K = 2Li ± 2L2 ± L3 - 2L4 - y ± 2L6 ± Lg . (66) 



At tree-level, Weinberg [75] determined that m.^ai/2 = 0.137 and 771^03/2 = —0.0687. 
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FIG. 8: The effective tt+K"*" scattering length times the reduced mass, Ht^k CL-iz+K+i't) as a function 
of time slice arising from smeared sinks. The solid black lines and shaded regions are fits with 1-a 
errors. The dashed lines on the tjItj ~ 350 MeV ensemble plot are an estimate of the systematic 
error due to fitting. 



The I = \ and / = i scattering lengths are related to crossing even and odd amplitudes by 



ai/2 = a"*" + 2a 
«3/2 = - a" 



0"iT+K+ 



(67) 



It is convenient to define the function F via a subtraction of the tree-level and one-loop 
contributions in order to isolate the counterterms, 



F 



2 [^^^TK a-K+K+\ + J- + X 2 — X ' 



(68) 



35 




FIG. 9: iIt^k <in+K+ vs. h-kk/ fn- The data points are the results of this lattice calculation, while 
the curve is the theoretical prediction at tree level in chiral perturbation theory [79]. The dark 
error bar is statistical, while the lighter error bar corresponds to the systematic error. The vertical 
dashed line denotes the physical pion and kaon masses. 



and at NLO this becomes 



phys\ 



- 2 ^ L^Kifty^) 



(69) 
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The dependence of F on m^r and rriK determines L5 and L^^x and, in turn, allows an extrac 
tion of 03/2 and ai/2- The numerical values of F and their errors are plotted in Fig, 
fitting a straight line to the values of F as a function of mk/m.,^ the counterterms 
Lt^k (renormalized at f^^^^) can be determined. 



By 
and 



As described in Section [IV A 2 the results of the lattice calculations are pruned, to give 
'fit A" and "fit B" . With this limited data set it is not practical to fit to the NNLO expres- 



TABLE V: Results of the NLO fits. The values of 171^(13/2 and correspond to their extrap- 

olated values at the physical point. 

FIT L5 X 10^ L^K X 10^ m^a3/2 rn^ai/2 xVdof 

A 3.83 ±0.49 3.55 ±0.20 -0.0607 ±0.0025 0.1631 ±0.0062 0.17 
B 2.94 ±0.07 3.27 ±0.02 -0.0620 ±0.0004 0.1585 ±0.0011 0.001 
C 5.65 ±0.02+°;^^ 4.24 ±0.17 -0.0567 ±0.0017 0.1731 ±0.0017 0.84 
D 5.65 ±0.02+°;^^ 4.16 ±0.18 -0.0574 ±0.0016 0.1725 ±0.0017 0.90 



sion [130] for the scattering length. However, it is important to estimate the uncertainty in 
the values of the scattering lengths extrapolated to the physical point that is introduced by 
the truncation of the chiral expansion at NLO. In our work on fx/ fn |136j a value of L5 was 
extracted as it is the only NLO counterterm that contributes. The numerical value obtained 
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FIG. 10: r vs. mx/mn.,^. The dark error bar on the data points is statistical, while the lighter error 
bar corresponds to the systematic error. The lines correspond to the four linear fits (A,B,C,D). 
The bars on the y axis represent the 1-a errors in the determinations of L5 = V{mK /m-K = 0) (see 
Ref. [133] ). (At 95% confidence level, these determinations are in agreement.) 



is only perturbatively close to its true value, as it is contaminated by higher-order contribu- 



tions. Therefore, by fixing the L5 that appears in Eq. (69) to the value of L5 extracted from 
Jk/ f-Ki an estimate of the uncertainty in both Lt^k and in the extrapolated values of the 
scattering lengths due to the truncation of the chiral expansion can be made. Specifically, 
L5 was sampled from a Gaussian distribution for a range of fx/ f-w values |136j and then 
Lt^k was fit. This fit is denoted "fit C", and the same fit but with the ~ 590 MeV data 
pruned is denoted "fit D" . The results of the four fits are given in Table |V] and plotted in 
Fig. [10} These lead to 

L^K = 4.16±0.18+°:26 , (70) 
and a prediction of the scattering lengths extrapolated to the physical point of 

ay2 = -0.0574 ± 0.0016+°:°°^4 
m^ai/2 = 0.1725 ± 0.00171°:°?! . (71) 

The central values and statistical errors were taken from fit D and the systematic error due 
to truncation of the chiral expansion was set by taking the range of the various quantities 
allowed by the four fits, including statistical and systematic errors. Fig. 11 shows the 95% 



confidence-level error ellipses associated with the four fits in the ai/2-'"^7r 0,3/2 plane. For 
purposes of comparison the current-algebra point [79j is included on the plot as well as 1-a 
error ellipses from analyses based on fitting experimental data using X"PT at NLO |127j 
and using Roy-Steiner equations |123] . As 1-a error ellipses correspond to 39% confidence 
level, one should be careful in finding discrepancy between the various determinations of the 
scattering lengths. The crossing-odd scattering length is of special interest as its corrections 
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FIG. 11: Error ellipses for the four fits (A,B,C,D) at 95% confidence level. The star corresponds 
to the current- algebra predictions (x-PT p^) from Ref. j79j . The 1-a error ellipses from a X"PT 
analysis at NLO [127 (denoted p'^) and from a fit using the Roy-Steiner equations [T1E\ are 

also shown. 



are protected by SU(2) chiral symmetry and are therefore of order and expected to be 
small [1371 [EH] 

Given how well the lattice data fit the NLO continuum X"PT formulas, it would seem 
that the 0(5^) discretization errors are comparable or smaller than the systematic error due 
to omitted 0{m^) effects in the chiral expansion. However, one should keep in mind that our 
determinations of, for instance, the low-energy constants L5 and L^^k are subject to 0(6^) 
shifts. In contrast with the 7r+7r+ and K~^K~^ scattering lengths, the mixed-action quantity 
Amix makes an exphcit contribution to the K^n^ scattering length [5 mil39j . While this adds 
an additional unknown contribution to this process, a MAx-PT analysis of nK scattering, 
including lattice data from the fine MILC lattices {b ~ 0.09 fm), as well as the determination 
of Amix [ 140J , will be able to address this source of systematic error quantitatively. 

We anticipate that with improved statistics, together with calculations on lattices with 
smaller lattice spacings, the theoretically-predicted regions for ^,^03/2 and m.,^ai/2 can be 
further reduced beyond those shown in Fig. 11 These regions can then be compared with the 
expected measurements from K~^tt~ atoms, to provide an exciting test of hadronic theory. 

We also note that (1) there are calculations underway to calculate the scattering lengths 
in both isospin channels directly with quenched QCD |141] : (2) an indirect method which 
uses scalar form factors for exclusive semileptonic decays, calculated with nj = 2 lattice 
data, gives ai/2 = 0.179 ± 0.017 [11211113] . 
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We thank Heiri Leutwyler for emphasizing this point to us. 
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C. KK Scattering 



1. Introduction 

Strange hadrons may play a crucial role in the properties and evolution of nuclear material 
under extreme conditions [1] . The interior of neutron stars provide one such environment in 
which the densities are high enough that it may be energetically favorable to have strange 
baryons present in significant quantities, depending upon their interactions with non-strange 
hadrons. Further, it may be the case that a kaon condensate forms due, in part, to strong 
interactions between kaons and nucleons [144]. Unfortunately, the theoretical analysis of 
both scenarios is somewhat plagued by the limited knowledge of the interactions of strange 
hadrons with themselves and with non-strange hadrons. 

Heavy-ion collisions, such as those at the BNL Relativistic Heavy Ion Collider (RHIC), 
also produce nuclear material in an extreme condition. Recent observations suggesting the 
formation of a low- viscosity fluid are quite exciting as they provide a first glimpse of matter 
not seen previously. The late-time evolution of such a collision requires an understanding 
of the interaction between many species of hadrons, not just those of the initial state, 
including the interactions between strange mesons and baryons. While pion interferometry 
in heavy-ion collisions is a well-established tool for studying the collision region (for recent 
theoretical progress, see Refs. jl45[ll46]ll47j ). the STAR collaboration has recently published 
the first observation of neutral kaon (K^) interferometry [148J . In the analysis of K^-K^ 
interferometry, the non-resonant contributions to the final state interactions between the 
kaons were estimated using three-flavor X"PT, the low-energy effective field theory of QCD. 
Given the sometimes poor convergence of SU (3) l ® SU (3)^ x-PT due to the relatively large 
kaon mass compared to the scale of chiral symmetry breaking (A^ ~ 1 GeV), particularly in 
the baryon sector, it is important to be able to verify that the non-resonant contributions 
to i^'ii'-scattering are indeed small, as estimated in X"PT. To date, there have been no 
experimental determinations of the 1 = 1 KK scattering length, but recently it has 

been calculated at NLO in x-PT |31i| . 

2. Mixed-Action x-PT One Loop 

In Ref. [31], the expression for the 1 = 1 KK scattering length was determined to NLO in 
X-PT, including corrections due to mixed-action lattice artifacts. As with the / = 2 tttt scat- 
tering length [SU], it was demonstrated that when the mixed-action extrapolation formula is 
expressed in terms of the lattice-physical parameters there are no lattice-spacing-dependent 
counterterms at 0{b'^), 0{h'^m\) or (9(6^). There are finite lattice-spacing-dependent cor- 
rections, proportional to 6^Ai, and therefore entirely determined to this order in MA^-PT. 
Again, as with the / = 2 tttt system, the NLO MA formula for itlkOj^kk does not depend 
upon the mixed valence-sea meson masses, and therefore does not require knowledge of the 
mixed-meson masses |140] . This allows for a precise determination of the predicted MA 
corrections to the scattering length. At NLO in MA^-PT, the scattering length takes the 
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FIG. 12: The effective K~^K^ scattering lengtli times tlie effective mj^+ as a function of time slice. 
The solid black lines and shaded regions are fits with 1-a statistical uncertainties. The dashed lines 
correspond to the statistical and systematic uncertainties added in quadrature. (See Ref. |l.l49j for 
data tables.) 
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(72) 



where the various coefficients, along with fh\ and mf^, can be found in Appendix 
E of Ref. [31]. To account for the predicted MA corrections, one can use Eq. (72) to 
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FIG. 13: mx+aK+K+ versus rn^+ / fj^+. The points with error-bars are the results of this lattice 
calculation (not extrapolated to the continuum) on both the coarse and fine MILC lattices. The 
solid curve corresponds to the tree-level prediction of X"PT, and the point denoted by a star 
and its associated uncertainty is the value extrapolated to the physical meson masses and to the 
continuum. 



directly fit the results of the lattice calculation. Further sources of systematic errors can 
be identified; higher-order effects in the chiral expansion, ^NNLo{'n^KO'^KK)'^ exponentially- 
suppressed finite- volume effects, ^Fvi.'^nKO'^KK)'^ residual chiral symmetry breaking effects 
from the domain-wall action, Am^^^(mxa/fii'); and the error in truncating the effective-range 
expansion with the inverse scattering length, ^rangeij^KCi^KK) ■ -^^^ ^ detailed discussion of 
these various sources of uncertainty, we refer the reader to Ref. [149]. 



3. Extrapolation to the Physical Point 



The results of the lattice calculations were pruned (defined in Section IV A2) to explore 
the convergence of the chiral expansion, giving fits A, B, and C 

fx) 



The results of these fits are collected in Table VI 



nm \m\ for the low- 
The 



energy constant 

extracted values of L^j^^ from each of the fits are consistent with each other within the 
uncertainties. In analogy with the comparison convention employed for '7r"*"7r+, the results of 
the lattice calculation are extrapolated to the physical values of / /x+ = 0.8731 ±0.0096, 
iT^K+ 1 fK+ = 3.088 ± 0.018 and rrir,/ fK+ = 3.425 ± 0.0019 assuming isospin symmetry, and 
the absence of electromagnetism. Taking the range of values of L^j^^ spanned by these fits, 
gives 



mK+aK+K+ = -0.352 ±0.016 , 32(47r)2Ljf^(/i = /k) = 7.1 ± 0.7 ^ 



(73) 
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TABLE VI: The results of fitting three-flavor MAx-PT at NLO to the computed scattering lengths, 
as described in the text. The values of mx+aK+K+ a-re those extrapolated to the physical (isospin- 
symmetric) meson masses and to the continuum. 



FIT 


32(47r)2L^=i(/;,) 


nif^+ax+K+ (extrapolated) 


xVdof 


A 


7.3(1)(4) 


-0.347 ±0.003 ±0.009 


0.22 


B 


7.3(2)(5) 


-0.347 ±0.004 ±0.011 


0.32 


C 


6.9(2)(6) 


-0.355 ±0.005 ±0.013 


0.14 



where the statistical and systematic errors have been combined in quadrature. The results 



are shown in Fig. 13 



D. Thoughts on Meson- Meson Scattering 

The lattice results for meson-meson scattering pose an interesting puzzle. In Fig. [6] one sees 
that the / = 2 tttt scattering length tracks the current algebra result up to pion masses that 
are expected to be at the edge of the chiral regime in the two-flavor sector. While in the two 
flavor theory one expects fairly good convergence of the chiral expansion and, moreover, one 
expects that the effective expansion parameter is small in the channel with maximal isospin, 
the lattice calculation clearly imply a cancellation between chiral logs and counterterms 



(evaluated at a given scale). As one sees in Fig. 13, the same phenomenon occurs in K^K^ 
where the chiral expansion is governed by the strange quark mass and is therefore expected 
to be more slowly converging. The tt^K^ scattering length (see Fig. |9| indicates similar 
behavior but one should keep in mind that this data is expected to contain larger lattice 
spacing corrections. This mysterious cancellation between chiral logs and counterterms for 
the meson-meson scattering lengths begs for an explanation. 



E. NN Scattering 

One of the ultimate goals of nuclear physics is to compute the properties and interactions 
of nuclei directly from QCD. The first, pioneering study of NN scattering with Lattice 
QCD was performed more than a decade ago by Fukugita et al [99j. This calculation was 
quenched and at relatively large pion masses, m^r^ 550 MeV. A fully dynamical Lattice 
QCD calculation of the NN scattering lengths in both the ^^o-channel and ^5*1 — '^Di-coupled- 
channels at pion masses of m.,^ ~ 350 MeV, 490 MeV and 590 MeV in the isospin-limit was 
performed within the last few years |150] . The dependence of the NN scattering lengths 
upon the light-quark masses has been determined to various non-trivial orders in the EFT 
expansion |151t 11521 1153] , and is estimated to be valid up to m-n ~ 350 MeV. The pion 
mass dependence of nuclear observables may also be explored using a lattice formulation 
of the low-energy EFT of nucleons [154J. The results of the Lattice QCD calculation at 
the lightest pion mass and the experimentally-determined scattering lengths at the physical 
value of the pion mass are used to constrain the chiral dependence of the scattering lengths 
from m^r ~ 350 MeV down to the chiral limit. 

For an arbitrary nucleus (or bound and continuum nucleons), of atomic number A and 
charge Z, there are naively {A + Z)\ {2 A — Z)\ contractions that must be performed to 
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produce the nuclear correlation function Therefore, in the ^5*0 channel there are 48 
contractions, while in the ^5*1 — ^-Di coupled channels system there are 36. The two-nucleon 
correlator that projects onto the s-wave state in the continuum limit is 

(±\ _ yijkl 
^NNV') — ^a/3ap 

5^(iVf(t,x)iV/(t,y)iVr(0,0)iV,''t(0,0)) , (74) 

where / denotes the isospin of the NN system and 5* denotes its spin, a,P,a,p are isospin- 
indices and i,j,k,l are Dirac-indices. The tensor X^^^l^ has elements that produce the 
correct spin-isospin quantum numbers of two-nucleons in an s-wave. The summation over 
X (and y) corresponds to summing over all the spatial lattice sites, thereby projecting 
onto the momentum p = state. The interpolating field for the proton is Pi{t,x.) = 
eabcUiit,x.) (^''^(t, x)C75(i^(t, x)), and similarly for the neutron. 

At the pion masses used in these calculations the NN scattering lengths are found to 
be of natural size in both channels, and are much smaller than the L ~ 2.5 fm lattice 
spatial extent. It is noteworthy that the scattering lengths at the heaviest pion mass are 
not inconsistent with the hghtest-mass quenched values of Ref. [99]. However, one should 
keep in mind the effects of quenching on the infrared properties of the theory |156] . 

A pion mass of m^r ~ 350 MeV is expected to be at the upper limit of where the EFT 
describing NN interactions will be valid |157[ 11581 1159[ I160[ 11611 1162j . The chiral extrapo- 
lation is performed with BBSvK power-counting |162] (=KSW power-counting |16Ul I161j ) 
and W power-counting |157l 11581 1159] in the ^S'o-channel and BBSvK power-counting in the 
^5*1 —^Di coupled channels. The Lattice QCD determinations of the light-quark axial-matrix 
element in the nucleon by LHPC [63] and its physical value are used to constrain the chiral 
expansion of qa- The lattice calculations of the nucleon mass and pion decay constant [78] 
— as well as their physical values — are used to constrain their respective chiral expansions. 
In addition to the quark-mass dependence these three quantities contribute to the NN sys- 
tems, there is dependence on the quark masses at next-to-leading order (NLO) from pion 
exchange, and from local four-nucleon operators that involve a single insertion of the light- 
quark mass matrix, described by the "1^2" coefficients [151^ 11521 1153j . The results of this 
Lattice QCD calculation constrain the range of allowed values for the D2^s, and consequently 
the scattering lengths in the region between m^r ~ 350 MeV and the chiral limit, as shown 



in Fig. 14 and Fig. 15 With only one lattice point at the edge of the regime of applicability 
of the EFT, a prediction for the scattering lengths at the physical pion mass is not possible: 
the experimental values of the scattering lengths are still required for an extrapolation to 
the chiral limit and naive dimensional analysis (NDA) is still required to select only those 
operator coefficients that are consistent with perturbation theory. The regions plotted in 
the figures correspond to values of Cq - the coefficient of the leading-order quark-mass in- 
dependent local operator - and D2 that fit the lattice datum and the physical value, and 
are consistent with NDA. In both channels the lightest lattice datum constrains the chiral 
extrapolation to two distinct bands which are sensitive to both the quark mass dependence 
of Qa and the sign of the D2 coefficient. As the lattice point used to constrain the EFT is at 
the upper limits of applicability of the EFT, we expect non-negligible corrections to these 



regions from higher orders in the EFT expansion. It is clear from Fig. 14 and Fig. 15 that 
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The sum over all contractions can be shown |155| to scale as A^, rather than ^ {Al^ 
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even a qualitative understanding of the chiral limit will require lattice calculations at lighter 
quark masses. 




FIG. 14: Allowed regions for the scattering length in the 5*0 channel as a function of the pion 
mass. The experimental value of the scattering length and NDA have been used to constrain the 
extrapolation in both BBSvK [l60l HSU [l62] and W [Wl [Ml [EH] power-countings at NLO. 
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FIG. 15: Allowed regions for the scattering length in the "^Si —^Di coupled-channels as a function 
of the pion mass. The experimental value of the scattering length and NDA have been used to 
constrain the extrapolation in BBSvK |162j power-counting at NLO. 

Without the resources to perform similar Lattice QCD calculations in different volumes, 
and observing that most energy-splitting are positive, it is assumed that scattering states 
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are found in each case. Calculations in larger volumes will be done in the future to verify the 
expected power-law dependence upon volume that scattering states exhibit. In addition to 
discriminating between bound and continuum states, calculations in a larger volume would 
reduce the energy of the lowest-lying continuum lattice states, and thus reduce the uncer- 
tainty in the scattering length due to truncation of the ERE. Further improvement would 
result from measuring the energy of the first excited state on the lattice, either with a single 
source or by using the Liischer- Wolff method. The lattice spacing effects in the present 
calculation appear at ~ (or exponentially suppressed 0{b) effects), and are expected 

to be small. However, the finite lattice spacing effects should be determined by performing 
the same calculation on lattices with a finer lattice spacing. A theoretical investigation of 
mixed-action EFT for NN scattering which would allow a continuum extrapolation remains 
to be carried through In addition to more precise Lattice QCD calculations through 
an increase in computing resources, formal developments are also required. In order to 
have a more precise chiral extrapolation, calculations in the various relevant EFTs must 
be performed beyond NLO. Furthermore, it is clear that lattice calculations at lower pion 
masses are essential for the extrapolation to the chiral limit, and will ultimately allow for a 
"prediction" of the physical scattering lengths. 

F. YN Scattering 

1. Introduction 

In high-density baryonic systems, the large value of the Fermi energy may make it en- 
ergetically advantageous for some of the nucleons to become hyperons, with the increase 
in rest mass being more than compensated for by the decrease in Fermi energy. This 
is speculated to occur in the interior of neutron stars, but a quantitative understanding 
of this phenomenon depends on knowledge of the interactions among the hadrons in the 
medium. In contrast to NN interactions, where the wealth of experimental data has allowed 
for the construction of high-precision potentials, the hyperon-nucleon (YN) interactions are 
only very-approximately known. Experimental information about the YN interaction comes 
mainly from the study of hypernuclei |164[ I165j . the analysis of associated A-kaon and E- 
kaon production in NN collisions near threshold |166t 11671 1168[ I169[ 11701 1171j . and hadronic 
atoms |172] . There are a total of 35 cross-sections measurements |173] of the processes 
Ap — > Ap, S^p — > An, T,^p — ^ S+p, T,^p S^p and S^p — > S°n, and unsurprisingly, 
the extracted scattering parameters are highly model dependent. The theoretical study of 
YN interactions is hindered by the lack of experimental guidance. The "realistic" potentials 
developed by the Nijmegen [1741 1175] and Jiilich |176l 11771 1178] groups are just two examples 
of phenomenological models based on meson exchange. These are soft-core potentials with 
one-boson exchange models of the NN interaction. Since SU{3) flavor symmetry is bro- 
ken by the differences in the quark masses, the corresponding couplings are not completely 
determined by the NN interaction and are instead obtained by a fit to the available data. 
In Ref. |174[ I175j . for example, six different models are constructed, each describing the 
available YN cross-section data equally well, but predicting different values for the phase 



See Ref. |34j for preliminary work in this direction. The masses and magnetic moments of the octet 
baryons have been calculated at NLO in MA^-PT in Ref. [163] . 
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shifts. The effective field theory approach jl79[ I180[ 11811 11821 1183] is less developed and 
suffers from a large number of couplings that need to be fit to the data. 

In view of the large uncertainties in the YN scattering amplitudes and their importance 
for modeling neutron stars and the study of hypernuclei, a first-principles QCD calculation 
of YN scattering is highly desirable. The only way to achieve this is through numerical 
calculations using Lattice QCD. Some of the present authors were part of a work |182j that 
outlined a program to address this issue with a combination of lattice calculations and the 
use of effective field theories. 

We have computed the low-energy s-wave phase shifts for YN scattering |184j in the ^5*0 
channel and ^5*1 —^Di coupled-channels at particular energies, using Liischer's finite-volume 
method [T6| HU |T9] . No attempt was made to extrapolate to the physical pion mass as it is 
likely that all but one of the data points lies outside the regime of applicability of the YN 
EFTs. 

An interpolating field of the form n x E~ (a simple product of the single baryon inter- 
polating fields) was used to determine the energy-eigenvalues of the s-wave strangeness = 1, 
isospin = I eigenstates in both spin channels, and an interpolating field n x A to determine 
the energy-eigenvalues of the s-wave strangeness = 1, isospin = ^ eigenstates in both spin 
channels. 



2. Results 



A typical example of the quality of the output of the lattice calculation can be seen in the 
effective mass plots for nA in the ^^o-channel, as shown in Fig. 16 The plateaus in the 
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FIG. 16: The effective mass plots for nA in the ^S'o-channel at pion masses of ~ 350 MeV (left 
panel), ~ 490 MeV (center panel) and m,r ~ 590 MeV (right panel). The straight line and 
shaded region correspond to the extracted energy shift and associated uncertainty. The dashed 
lines correspond to the statistical and systematic errors added linearly. 

effective energy plots persist for only a small number of time-slices. At small t there is the 
usual contamination from excited states whereas at larger t the signal-to-noise ratio degrades 
exponentially with t. The Dirichlet boundary at t = 22 introduces a systematic uncertainty 
due to backward propagating states. However, in practice, the statistical noise becomes a 
limiting factor at far earlier time slices and the boundary at t = 22 is not an issue for this 
calculation. We obtained a non-zero energy shift larger than the statistical error in ten of 
the correlation functions. 

It is not clear that we have been able to identify the ground states in all of the correlation 
functions, e.g. nS~ in the ^5'o-channel at ~ 490 MeV, and nA in the ^S'o-channel at 
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~ 490 MeV, as the statistics are not sufficient to determine whether the large-time 
behavior we observe is due to noise or due to the presence of any states with lower energy 
than those found. Indeed, it would be very exciting if there were states with lower energy, as 
they would likely be bound states (based on naturalness arguments and the exact Liischer 
relation) . 
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FIG. 17: Comparison of the lowest-pion-mass lattice results in each channel with a recently de- 
veloped YN EFT |183] (squares), and several potential models: Nijmegen [174] (triangles) and 
Jiilich |178] (diamonds). The dark error bars on the lattice data are statistical and the light error 
bars are statistical and systematic errors added in quadrature. 



3. Discussion 

The lightest pion mass at which a signal was extracted is at the upper limits of the regime 
of applicability of the effective field theories that have been constructed, thus precluding a 
chiral extrapolation. However, this work does provide new rigorous theoretical constraints 



on effective field theory, and potential model constructions of YN interactions. In Fig. [17 
we compare the lattice values of the phase shifts to recent EFT results jl83j (squares), and 
to several potential models: Nijmegen jl74] (triangles) and Jiilich jl78] . At face value these 
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results appear quite discrepant, however one should keep in mind that extrapolation to the 
physical pion mass will seriously alter individual contributions to the YN interaction. 

While the measurements of the momenta and phase shifts are unambiguous, their physical 
interpretation is not entirely resolved. Each of the phase shifts at the lowest pion masses are 
negative. Assuming that the observed state is the ground state in the lattice volume, this 
implies that the interactions are all repulsive. The raS" interaction in the ^5*1 —^Di coupled 
channels is strongly repulsive at ttItt ~ 490 MeV, while the interaction in the ^So-channel is 
only mildly repulsive. The opposite is found to be true for the nK systems at ~ 350 MeV, 
where the interaction in the ^S'o-channel is found to be strongly repulsive, while that in 
the ^5*1 — ^-Di coupled channels is mildly repulsive. However, there may be channels for 
which there exist states of lower, negative energies, some of which may correspond to bound 
states in the continuum limit. If such states are present, then we would conclude that 
the interaction is attractive, and that the positive-shifted energy state we have identified 
corresponds to the first continuum level. Current statistics are sufficiently poor that nothing 
definitive can be said about the existence of such states. 



G. Exploratory Quenched Calculations 

1. BB Potentials 

Energy-independent potentials can be rigorously defined and calculated for systems com- 
posed of two (or more) hadrons containing a heavy quark in the heavy-quark limit, uiq — > oo. 
This is interesting for more than academic reasons as the light degrees of freedom (dof) in 
the B- meson have the same quantum numbers as the nucleon, isospin-^ and spin-|. As 
such, the EFT describing the interactions between two B-mesons has the same form as 
that describing the interactions between two nucleons, but the counterterms that enter 
into each EFT are different. Therefore, a deeper understanding of the EFT description 
of nuclear physics can be gained by Lattice QCD calculations of the potentials between 
B-mesons. A number of quenched explorations of these potentials have been performed 
previously [UHl 1861 [Wl [Ml [1891 [M [1911 M [193] , but there was little evidence for any 
potential simply due to limited computational resources. We returned to this problem |194j 
and computed the potential between two B-mesons in the four possible spin-isospin channels 
neg lecting - mixing) in relatively small volume DBW2 lattices with L ~ 1.6 fm, 
with a pion mass of ~ 403 MeV, and lattice-spacing of 6 ~ 0.1 fm. The calculation was 
quenched and the naive Wilson action was used for the quarks. At this relatively fine lattice 
spacing, much finer than previous calculations, were able to extract a non-zero potential, 
but the small volume meant that the contributions to the potential from image B-mesons 
(periodic BC's) were visible. 

Constructing the t-channel potentials, defined via the quantum numbers of the exchange 
particles, in keeping with nuclear physics tradition, isolated statistical fiuctuations into the 
channel associated with the "a" -meson, leaving the channels with the quantum numbers of 
the TT, p and ui with relatively small statistical errors. The potentials are shown in Fig 



Given the uncertainties in the potentials, and the number of counterterms that appear 
in the EFT describing the long- and medium- distance interactions between the B-mesons, 
it was possible to make only a parameterization of each potential beyond the leading light- 
meson contribution. Since only the longest range contribution to the potential in each 
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FIG. 18: The t-channel potentials between two B-mesons in the finite lattice volume. 
Va^^ , Vr^^ , Var '^ , and V^^^ coriespond to the potentials in the exchange-channels with spin-isospin 
of (J,/) = (1,0), (0,1), (1,1) and (0,0). 




FIG. 19: Fits to the finite-volume isovector t-channel potentials between two B-mesons. The 
dashed lines correspond to the finite-volume fits to the lattice data, and the solid curves are the 
infinite- volume extrapolations. 

channel can be identified, we fit our results at large separations, |r| > A^^, using the finite- 
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volume versions of the simplified infinite-volume potentials, 

V^^\y) ^ " ^ + a' , (75) 

'^^ ^ ' 247r/2 |r| ^ |r| ' ^ ^ 

, , lrl->no p-"^p|r| p-^xkl 

V;(-)(r) + . (76) 

^ ^ ' An |r| |r| 

The short distance forms in the above equations are entirely model dependent and are the 
simplest forms that we could find that provide a reasonable description of the data. Using 
the measured values and uncertainties of m.^ and nip and the physical value of f-^ we first 
determine the couplings g and gp by setting = a'^ = and fitting the finite-volume 
potentials at the two largest separations.^^ These fits are shown by the dashed red curves 



in Fig. [19] and the resulting couplings are found to be 

= 2.17 ±0.08, c/ = 0.57 ±0.06. (77) 

Having determined these parameters, we reconstruct the infinite-volume potentials that are 
shown in the figure as the solid red lines. 

Clearly the lattice calculations that exist of the potentials between B-mesons must be 
viewed as nothing more than exploratory. Further, the analysis of the output of the cal- 
culations should be model-independent in order to impact our understanding of the NN 
interaction. 

2. J/ip-Hadron Scattering 

The interactions between quarkonia and the light hadrons is interesting from both the theo- 
retical and experimental standpoints. Theoretically, near the heavy-quark limit, the quarko- 
nia are compact objects for which a multipole expansion of the chromo-electric and magnetic 
fields can be performed, with the leading interactions scaling as f'^^g, where Tqq is the ra- 
dius of the QQ state. This has lead to predictions for the binding of such states to nuclear 
matter [195j and interactions with light nuclei |196j . On the experimental side, the number 
and distribution of quarkonia observed in heavy ion collisions as a function of nuclear size 
provides important information on both the production and attenuation of such states in nu- 
clear matter under extreme conditions. The first quenched lattice studies of the low-energy 
interactions of charmonium with a tt, p or nucleon have been recently performed |197j . This 
work is quite encouraging, and we look forward to fully-dynamical calculations at smaller 
pion masses. 

V. N-BODY INTERACTIONS 

The simplest multi-hadron systems consist of n pseudoscalar mesons of maximal isospin. 
Although such systems are of limited phenomenological interest for n > 3 (three pion in- 
terferometry is currently a topic of interest in heavy- ion collisions |198[ I199[ 1200] ). they 
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Simple fits using the infinite- volume long range behavior were considered in Ref. |185j . 
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FIG. 20: The effective energy plots associated witli tlie correlation functions of n 
with the pion rest mass removed for a pion mass of ~ 290 MeV. 



2,3, ..,12 vr+'s 



serve as a testing ground for more complicated many-body systems 
first Lattice QCD calculation of multi-hadron {n > 2) interactions 



We performed the 
by calculating 



the ground-state energies of tt+tt^ 



7r+7r+7r+7r+ up to (tt^ 



\12 



in a spatial volume of 



V ~ (2.5 fm)^ with periodic boundary conditions. These systems serve as an ideal laboratory 
for investigating multi-particle interactions as chiral symmetry guarantees relatively weak 
interactions among pions, and multiple pion correlation functions computed with Lattice 



QCD do not suffer from signal to noise issues (as can be seen in Fig. 20) that are expected 



to plague analogous calculations in multi-baryon systems, as discussed in section |IIE[ The 
TT+TT"*" scattering length is extracted from the n-pion {n > 2) systems with precision that is 
comparable to (and in some cases better than) the n = 2 determination |115j . 

Using the expression for the energy shift of ?T,-7r"'"'s in a finite volume due to their two- 



and three-body interactions, as given in Eq. (14), the three- vr"*" interaction 



=L -L 



1 



6 ( — ) X 

.71 L, 



(78) 



L . 



is determined at four different quark masses. The extracted values of r]^ in units of the NDA 
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FIG. 21: The 7r+7r+7r+ interaction, 773, in units of the NDA estimate of 1/(7717^ /^) as a function 
of m^//^. 



estimate of 1/(^7^ /^) are shown in Fig. 21 as a function of niT^/ Jt^. There is clear evidence 
for a three-TT"*" interaction, and further, it is of a magnitude that is consistent with NDA. 

The study of systems with arbitrary numbers of vr+'s in a fixed volume enables the study 
of the isospin chemical potential as a function of isospin density. For the systems containing 
less than twelve tt+'s with the coarse MILC lattices with L ~ 2.5 fm, the results are shown 



in Fig. 22 [23] . This is a first step toward a general study of multi-hadron systems. 



VI. RESOURCE REQUIREMENTS FOR FURTHER PROGRESS 

To understand more concretely the resources that will be required to compute the interac- 
tions between two or more nucleons, we determined the uncertainty in the nucleon-nucleon 
scattering length, at a scattering length of a = 2 fm (independent of the pion mass) as a 
function of the computational resources available to this program. Our estimates, shown in 



Fig. 23, are for computations on lattices with a lattice spacing of 6 ~ 0.125 fm and a spatial 
lattice extent of L ~ 5 fm. The following procedure was used to determine the uncertainties 
in the scattering length: 

1. The results obtained by NPLQCD indicate that the uncertainty in the scattering 
lengths depend (approximately) exponentially upon the pion mass. The argument 
and coefficient of the exponential were fit to these results. 

2. The computational requirements for propagator generation scale as Tcpu = {A + 
B/mg)V. A and B were determined by the timings on the b ~ 0.125 fm, 20^ x 64 
MILC lattices for domain-wall propagator generation. The lattice volume is in lattice 
units. 

3. The value of the energy splitting between the two nucleons in the volume and two 
isolated nucleons in the volume was tuned to produce a scattering length of 2 fm. 
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FIG. 22: The isospin chemical potential as a function of the number of pions (equivalent to 
the isospin density pj) at a fixed volume with and without the contribution from the vr^Tr+vr^- 
interaction, T/3 . The solid red (lighter) curves and bands result from the analytic expression for the 
energy of the ground state in the large volume expansion, Eq. (14), using the fit values for a^+^+ 
and r/g and their correlated uncertainties. The solid blue (darker) curves and bands are similarly 
the results for the fitted value of a,r+7r+ Vs—^- The dashed curve corresponds to the leading 
order prediction of X"PT. 



4. The statistical uncertainty on the projected energy-splitting was determined by the 
number of propagators that could be generated with the given computational resources. 
Synthetic data was then processed using the jackknife procedure with the exact Luscher 
formula to produce the uncertainty in the scattering length. 

5. We have not included the effect of the improved signal in going to larger volumes, 
possibly a 1/ ^/V factor, as we have no data to extrapolate from. 

The quoted resource requirements are for propagator generation ONLY, and do not include 
the resource requirements for lattice generation. Each data point requires the resource 
indicated in the legend of Fig. 123} Therefore, the computational resources required to 



generate the data indicated by the purple squares in Fig. [23] with the precision that is shown, 
is 0.12 Pfiop-yrs. We assume that the lattices will be generated by others in the USQCD 
collaboration, and will be made available (in a timely fashion). Further, these resources are 
for the generation of Domain- Wall fermions. If, instead. Clover- Wilson fermions are used on 
dynamical Clover- Wilson lattices, the resource requirements for propagator generation are 
reduced by a factor of ~ 10. Additional factors of ~ 5 — 8 (included in these estimates) are 
gained using the newly developed Incremental and EigCG deflation method \201\ . Provided 



that gauge fleld conflgurations become available, the task of computing NN scattering lengths 
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FIG. 23: The expected uncertainty in nucleon-nucleon scattering in the —^Di coupled channels 
for a 2 fm scattering length as a function of the pion mass for domain wall valence quarks on the 
coarse MILC lattices (6 ~ 0.125 fm) for given computational resources (there is a further reduction 
by an order of magnitude for clover quarks in the valence sector). Each data point requires the full 
computational resource, i.e. to generate the complete set of six triangular (green) data requires 
12 Tflop-yrs. 

and phase-shifts from lattice QCD could be accomplished with present-day computational 
resources. 



VII. CONCLUDING REMARKS 

Lattice QCD, when combined with effective field theory, is presently able to provide quan- 
titative information about the interactions among hadrons directly from QCD. During the 
last five years a number of two-hadron processes have been investigated in fully-dynamical 
lattice QCD, including both meson-meson and baryon-baryon processes. While, it is fair 
to say that the study of hadronic interactions with Lattice QCD is still in its infancy, the 
computational resources that will likely become available during the next five years will 
allow for a complete study such systems. This will amount to nothing short of a revolution 
in nuclear physics. 

The study of systems containing baryons is much more computationally intensive than 
that of systems containing only pions, due to the exponential degradation of the signal- 
to-noise in (multi-)baryon correlators. This is presently a serious roadblock for the direct 
study of nuclei and nuclear interactions with Lattice QCD. Further, the study of systems 
such as TT+vr" — > vr+vr" are presently out of reach due to the need to calculate disconnected 
diagrams, which require approximately an order of magnitude increase in computational 
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time over the connected contributions. 

For the calculations that are of interest to the nuclear physics community, it is important 
to step back and consider which lattice calculations are the most important to perform, and 
what the criterion "important" actually means! There may be some use in performing a high 
precision calculation, say ~ 0.1%-level, of a mesonic quantity, such as /^r, which is already 
known to high precision experimentally, but there may be equivalent use in performing a low 
precision calculation, say at the ~ 20%-level, of the weak pion-nucleon coupling constant, 
^i-AfiV' ^h^t expected to provide a significant contribution to parity violating observables in 
many light nuclei and for which there is still a substantial amount of uncertainty as regards 
its value. This is clearly something for future discussions and we do not attempt to resolve 
this issue in this review. 

The near future will see remarkable progress in this field. One can expect calculations 
in large volumes, small lattice spacings and at the physical quark masses, and also the 
inclusion of electromagnetism. With sufficient resources dedicated to this area of Lattice 
QCD, the interactions between two and three hadrons will be calculated at high precision, 
with complete control over all systematic uncertainties that arise. At that time, calculations 
of processes relevant to areas of nuclear physics, such as the interactions between hyperons 
and nucleons that impact the equation of state of nuclear matter in the interior of neutron 
stars, or the interactions between multiple nucleons, will be calculable directly from QCD. 
However, it is crucial to appreciate the fact that in order for this program in multi-baryon 
physics to be realized a significant amount of computational power will need to be directed 
exclusively at multi-baryon calculations. 

Finally, in direct analogy with the experimental programs, it is not sufficient to have 
just one calculation of any given quantity. Lattice calculations with different discretizations, 
lattice volumes, lattice spacing, performed by different lattice collaborations are required in 
order to have confidence in the "Lattice QCD" value of the quantity. This concerted effort 
does not yet exist for any scattering process, and we look forward to such a coherent effort 
taking shape in the near future. 
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